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Abstract. We consider the Hele-Shaw problem in a randomly perforated domain with zero Neumann 
boundary conditions. A homogenization limit is obtained as the characteristic scale of the domain goes 
to zero. Specifically, we prove that the solutions as well as their free boundaries converge uniformly 
to those corresponding to a homogeneous and anisotropic Hele-Shaw problem set in R d . The main 
challenge when deriving a limit lies in controlling the oscillations of the free boundary, this is overcome 
first by extending De Giorgi-Nash-Moser type estimates to perforated domains and second by proving 
the almost surely non-degenerate growth of the solution near its free boundary. 



1. Introduction 

For a stationary ergodic set O C M d and e set O e :— eO e we consider a nonnegative random function 

ii f :O e xl + -> K, 

which has compact support for t = and solves the free boundary problem 

-Au e = 1 in {u e > 0}, 

d n u e =0 on<90 £ , (1.1) 

u\ =|Vu e | 2 on c^-jV > 0}, 

where d n denotes the outer normal derivative on dO e and <9" denotes the boundary operator which acts 
on {it e > 0} relative to O e (oj), where to is the sample variable. The first two equations say u e solves an 
elliptic problem for fixed t and the third says {u e > 0} is expanding with normal velocity given by |Vu e |. 

This is the well known one-phase Hele-Shaw problem, posed in a random domain O e with Neumann 
boundary conditions. In this note, we shall prove a homogenization result which essentially says (see 



Theorem 2.6 for the exact statement) that almost surely u e and the free boundaries d u {u e > 0} converge 
uniformly to u and d{u > 0}, where solves the (deterministic) homogeneous problem 

f -div(AVu) =fi m{u>0}, 
1 u t =Q{Vu) on 9{ti > 0}. 



Here A, known as the effective diffusivity, is defined in terms of O (see Section 7, ( |7.3[ )). The constant 
/x represents the fraction of volume occupied by O and Q is the quadratic form Q(p) = (Ap ■ p)//i. Note 
that u e is defined in a different domain of M. d for each different e, but the putative limit u is defined 
everywhere (and in particular in O e for each e). Therefore, we define uniform convergence of u e to u as 

lim \\u(-,t) - u £ (-,i)|| L -(o ) = P-a.s. V t > 0, 
and similarly for the uniform convergence of the free boundaries (see Section [5] for precise definitions). 
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1.1. Motivation. The Hele-Shaw flow (for d = 2) describes movement of an incompressible, viscous 
fluid confined between two narrow plates, as first described in the seminal paper of Hele-Shaw |25j . In 
this context, u represents the pressure inside the cell, and our model arises as the surface tension zero 



limit when the plates trap a layer of some granular medium O e . Also for d = 2 one may see (1.1 1 as 
the evolution of a growing droplet which is spreading in a two dimensional labyrinth with walls lying 
vertically on top of dO e . The height of the droplet at time t is given by u e (.,t), thus the base of the 
droplet would be the "positive phase" {u e > 0}. Under this interpretation, the condition Au e = — 1 is a 
linearization of the constant mean curvature condition for the droplet, and the zero Neumann condition 
of u e just says the droplet meets the walls of the labyrinth at a right angle. Then, our main theorem 



(Theorem 2.6| states that in the limit e — > the droplets behave as if they were spreading on a free 
surface except their shape and motion law are changed, which typically become anisotropic (that is, the 
effective diffusivity A might not be diagonal, see Section [7]). 

Let us note however, that the main motivation for this work is not the justification of the physical 
validity of the homogenized limit for the droplet interpretation. Instead, the goal is to understand 
in a "simple" case, the regularity properties of diffusive free boundary problems when the underlying 
geometry is strongly irregular and random. By "strongly irregular" we mean that the randomness does 
not take the form of a diffusion operator in M. d with random uniformly elliptic coefficients, but instead 
that we are dealing with a random domain O e . A by product of our approach are uniform regularity 
estimates and pointwise bounds for solutions of linear elliptic equations in randomly perforated domains, 
which are new and of independent interest. 

Remark 1.1. A model which would be more physically realistic (for the situation of the droplet) is the 
one where we require —Au e = X(t) at each time, with A(i) > being picked so that the volume of the 
droplet is preserved over time. Our analysis could be extended to this case if we had a well-posedness 
theory for the problem for fixed e ( with some continuity estimate for X(t)), this is however not available 
at the moment. For more on droplet models, see [2TT [24l 12] and references therein. 



1.2. Background and contributions of this work. The evolution problem ( |1.1[ ), as well as the broad 
study of elliptic equations in O e is of interest both in the context of homogenization of partial differential 
equations as well as in percolation theory and analysis of (even discrete) stochastic growth models. First 
we review some of the relevant PDE and probabilistic ideas and how they relate to each other, followed 
by a discussion of our main results, which will be stated in detail in Section [2.3| 



Linear homogenization in random media. Recall that the random domain O e can be thought of as 
resulting from a degenerate limit of linear and uniformly elliptic operators in M d . Heuristically, for fixed 



e the evolution problem (1.1 ) is obtained as the limit as 5 — > of a family of Hele-Shaw problems posed 



in M. d with a random diffusivity matrix 



-dW(A s Vu s ) =1 in {u s > 0}, 

(ug) t = (AsVug^us) on 3{nj > 0}, 



where A$(x) is the matrix given by 

AAx) 



i if x e o e 

61 if x £ O e 



so in some sense, we are dealing with a degenerate elliptic problem in M. d where the diffusivity matrix 
field has the form A(x, u) = A(x,uA\q^(x) for some uniformly elliptic matrix A(x, oS). 

Let us focus for a moment on the theory of homogenization for (uniformly) elliptic operators in 
divergence form. Namely, (fixed) boundary value problems of the form 

f -div(A(f,w)Vu) =/ mDcR d , 
\ v = g on dD, 



Quasistatic Droplet Percolation 



3 



for a uniformly elliptic matrix process A(x,w), fixed domain D C R d and functions /, g. Then linear 
homogenization consists in finding whether the solutions to the above boundary value problems converge 
as e — > to the solution of a boundary value problem for a deterministic homogeneous operator. 

The first homogenization result was obtained by Tartar [3J Chapter 1, Section 3] in the case where 
A(x,ui) = A(x), for a periodic and uniformly elliptic A(x). Later, both Kozlov [34] and Papanicoloau- 
Varadhan .39J proved independently that homogenization holds for a stationary ergodic matrix process 
A(x,ui) (still assumed uniformly elliptic). Much later, the case where A is degenerate elliptic and 
periodic, that is A(x, lo) = Iq and O periodic (with R d \ O connected and with Lipschitz boundary) was 
studied by Cioranescu-Paulin |14j and Cioranescu-Donato |13j . The nonlinear theory in divergence form 
was done by Dal Maso - Modica [T3], see (SHUTS] for more complete references on all these topics. Most 



relevant to the analysis of (1.1) is the homogenization of the Neumann problem for the Laplacian in 
a stationary ergodic domain, the homogenization for this boundary value problem has been considered 
previously by Jikov |42j . 

The works of Kozlov and Papanicoloau-Varadhan already provide the almost sure convergence to 
the homogenized limit in L 2 , moreover, they show that the effective equation will be uniformly elliptic 
with the same ellipticity bounds as the heterogenous problems. In this case the matrix A(x, w) is 
uniformly elliptic, so one can use De Giorgi-Nash-Moser theory to get uniform interior Holder continuity 
of the solutions, guaranteeing the convergence is not just in L 2 but in fact locally uniform. In contrast, 
estimates of Dc Giorgi-Nash-Moser type are not readily availablein a random domain O t , and the question 
of (locally) uniform convergence has not been resolved yet in this case. 

Homogenization of free boundary problems. These problems take the previous theory as a starting 
point, as it provides at least the effective equation for the limiting free boundary problem away from 
the moving interface. The challenge is to identify the homogenized free boundary condition, which is 
not only a nonlinear condition but is also imposed on a lower dimensional (and often singular) set. A 



uniformly elliptic version of (1.1) was first achieved by Kim-Mellet in 28 j . They used the monotonicity of 
this one-phase problem as well as the uniqueness theory from [29] to show that the integral of a viscosity 
solution with respect to time solves a certain obstacle problem, something that was expected but had 
only been proved for classical solutions |18) . Using this correspondence, the homogenization of the 
original problem was reduced to the homogenization of an obstacle problem. Then, the homogenization 
of the obstacle problem can be handled since pointwise estimates (based on Harnack's inequality and 
barriers) control the free boundaries uniformly as e — > 0. 

There are also some results for uniformly elliptic one phase problems when one does not have any 
divergence structure (in which case the approach above fails), however, they are restricted to the periodic 
setting [3T]. These results do not rely on writing an auxiliary obstacle for the time integral, instead the 
approach is inspired by a different method introduced in [11] for stochastic homogenization of fully 
nonlinear elliptic equations. It is worth noting that this method also involves an auxiliary obstacle 
problem, although this is completely unrelated to the obstacle problem above. 

Related results in homogenization of free boundaries include singular limits for flame propagation 
problems, in particular recent work of Caffarelli, Lee and Mellet for the periodic case[TU] and the one- 
dimensional random case [3] ■ A time independent problem for droplets was considered by Mellet and 
Nolen{37], and homogenization of the one-phase Stefan problem was studied by the Kim and Mellet [32 . 

Remark 1.2. Another example of stochastic homogenization of a free boundary problem is the obstacle 
problem with a random highly oscillatory obstacle [5] . Although we also deal with an obstacle problem in 
the present work (as an auxiliary problem, see Section|4|, here we have a random domain with Neumann 
boundary conditions, a very different situation from the one treated in [5]. 

Percolation and growth models. From the probabilistic side, the Hele-Shaw problem in R d is the 
continuum counterpart (and often shown to be the continuous scaling limit) of several growth models 
taking place in the lattice Z d (such as the abelian sandpile) or internal diffusion limited aggregation 
(DLA) |35j . Keeping with this probabilistic interpretation of Hele-Shaw, our problem in a random 
domain is then a continuous analogue of internal DLA in a percolation cluster in TL d . 
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For a proper introduction to (discrete) percolation see Grimmet [23], see Barlow [3] for a summary 
of many of the results and analytical issues related to random walks in percolation clusters in Z d . The 
survey by Biskup [5] discusses many recent results and standing open questions. 

In the context of (continuous) random walks in infinite percolation clusters C C Z d , the questions cor- 
responding to homogenization are concerned with the long time behavior of the walk and the continuum 
limit when we let spatial scale of the lattice go to zero. For either case, obtaining pointwise bounds on 
the transition density of the random walk (i.e. the heat kernel for the percolation cluster) is essential, in 
other words, Nash-Aronson type bounds for the heat kernel are required. Equivalently, one is interested 
in proving De Giorgi-Nash-Moser type estimates in the infinite percolation cluster. 

However, it is well known that such pointwise estimates fail for supercritical percolation. In fact, in 
general, both continuum and discrete settings, a (uniform in cj) Sobolev inequality is not expected to 
hold. Heuristically, the infinite percolation cluster will always have some (maybe far away) region that 
is poorly connected, which could trap the random walk for a very long time, and this is reflected by 
the heat kernel. In geometric terms, the poor connectivity obstructs the validity of the isoperimetric 
inequality everywhere (and thus the Sobolev embedding and Poincare's inequality). 

This difficulty has been overcome, albeit only in the discrete case. The literature on such results is 
much too vast, we refer to the introduction in Barlow [3] as well as Biskup's survey [5] for a proper list 
of references. Also in [3J, a proof of non uniform bounds (in oj) is obtained by combining analytic and 
probabilistic ideas. Such bounds say that after waiting for some time (for a random time T x ) at a given 
point x <E % d the transition densities become Gaussian. The estimates by Barlow [3J go along the lines 
of the Fabes-Stroock proof of the Nash's inequalities j20j (for a deterministic equation). One may say 
that from all the approaches to De Giorgi-Nash-Moser theory, the one that has been most adaptable to 
the setting of discrete percolation is that of Nash. 

Remark 1.3. The validity of pointwise bounds on the Green's function (or heat kernel) for a domain 
O suggests that it enjoys some regularity properties. More precisely, if the Green's function of a given 
O is pointwise comparable with \x — y\ 2 ~ d , then not only does the Harnack inequality follows almost 
immediately, one can also prove a Sobolev inequality for functions in H (O). Compare this with the 
result that Gaussian bounds for the heat kernel (for a uniformly elliptic operator in divergence form) 
hold if and only if we have a certain doubling property of the volume and the Poincare inequality holds 
for all balls (see [221 [40] for Riemannian manifolds and [17 j for graphs, see also Appendix [B] ). 

This lack of uniform regularity estimates in C mirrors the lack of regularity estimates and uniform 
convergence for linear homogenization in stationary ergodic domains O e . This becomes a serious issue 



in the analysis of (1.1), since as hinted at earlier, one typically uses Harnack inequalities and regularity 
estimates to control the oscillations of the free boundary in the obstacle problem. If one is to carry 
out an approach similar to that of [28] . then one needs to develop an elliptic regularity theory in O e . 
We manage to do obtain uniform estimates in the continuum case under the extra assumption that O e 
is a randomly perforated domain (Assumption [2| , at the expense of not dealing with the continuum 
analogues of supercritical percolation models (such as the Boolean model). 

Contributions of this work. In this paper we consider randomly perforated domains O e , obtained by 
removing from R d a countable number of random isolated domains with Lipschitz boundaries. These 
go beyond the periodic domains studied in the linear homogenization literature, but do not cover many 
important percolation models (see the examples in Section[2]). On the other hand, we only need stationary 
ergodicity and do not make any assumptions about independence of the random variables. 



The main result for such domains O e concerns the free boundary problem (1.1), which we show 



homogenizes to (1.2), the convergence being almost surely uniform both for u e and the free boundaries 



9^{u e > 0} (see Theorem 2.6). This result also covers the one in [55], giving a simpler proof of the 
uniform convergence of the free boundaries. Moreover, our proof of homogenization has as a by product 
the homogenization of the usual obstacle problems in O e . 
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As a necessary step for the proof of this theorem we also prove uniform Holder estimates and Harnack 
inequality for harmonic functions in subsets of O e that have zero Neumann data on dO e (see Theorem 



2.7), such estimates, as well as the main homogenization theorem, are new even for the periodic case. 

Finally, we give a new (and shorter) proof of the homogenization result due to Jikov (see [32] also 
Chapter 8]) that is closer to the proofs in the uniformly elliptic case, at the expense of assuming that O e 
is a random perforated domain (see Assumption [2| . In particular, we use the regularity of the domain 



to characterize the space of correctors in O e (see Lemma 7.3). Moreover, we use the strong convergence 



granted by the uniform regularity estimates to prove that the effective diffusivity is strictly positive. 

Remark 1.4. It can be argued that the assumption of a perforated domain is rather strong, however, 
it is hoped that the machinery developed here combined with adequate isoperimetric inequalities for 
subsets of O e can be pushed to derive regularity bounds (Harnack, Nash-Aronson, etc) in continuum 
and discrete percolation using an approach closer to De Giorgi's (Nash's approach has been so far the 
most successful one in the discrete percolation literature). This will be explored in future work. On the 
other hand, as a payback for our stronger regularity assumptions, our elliptic estimates are uniform (in 
lu) and we apply them to get pointwise bounds for the Green's function in O e (see Appendix [Bj. 

1.3. Outline of the paper. The precise statement of the main homogenization theorem, as well as an 
in depth discussion of the strategy of the proof shall be deferred to later in Section [2] after we review the 
notation and basic definitions relevant to free boundary problems in random domains. 

The assumptions on O e and the main theorems are stated in Section [273] In Section [2~4| we discuss 
examples of random domains where they apply. In Section |2.5| we discuss the overall plan of the proof. 

Section [3] deals with the regularity of solutions to elliptic equations in O e with zero Neumann data 
on dO e , including Holder estimates and Harnack inequalities with bounds that are uniform in e. 

In Section [4] an important barrier is constructed, it is a non-negative function O e with quadratic 
growth at infinity (it plays the role that paraboloids play in M. d ) which has zero Neumann data on dO e . 
This construction is non-trivial and relies on the results of previous section. This barrier is used to 
show several important properties of the free boundaries d{u e > 0}. In particular, we show they cannot 
develop too many thin tentacles as e — > (i.e. in average they do not oscillate with large amplitude). 
In Sections[5]and[6]we combine the control on the free boundary oscillations with tools from stochastic 



homogenization and viscosity solutions to prove that (1.1) homogenizes. 

The new proof of homogenization of linear elliptic problems in random perforated domains is done in 
Section [7] Finally, in the appendices we prove the existence of weak solutions to the auxiliary obstacle 



problem derived from (1.1), and obtain uniform bounds for the Green's function in O e 



Acknowledgements. The authors would like to express their gratitude to Antoine Mellet for inter- 
esting discussions which prompted this research and to Marek Biskup for many helpful discussions on 
percolation. 

2. Free boundary problems in stationary ergodic domains 

2.1. Random perforated structures. We review some basic concepts and notation from the theory 
of stochastic homogenization for linear elliptic operators in divergence form [26) . 

Definition 2.1. An Ergodic Dynamical System is a pair where first we have a probability space (f2, J 7 , P) 
and second we have an ergodic action on ft by R d . The latter means we are given maps T x : — > O, 
x € M. d preserving the measure P, satisfying the group property T x+y = T x o T y , in a manner so that 

(x,u) -> T x u> 

is a measurable map from E'' x n (equipped with the product er-algebra) to f2, and such that the only 
events F 6 T which are preserved by all the maps T x are those with P(F) = or F(F) = 1. 

For the next definitions we take an ergodic dynamical system as given, we give some examples below. 
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Definition 2.2. Given a probability space (0, P) a real valued process is a measurable function 
/ : M. d x ft — > R, we will say / is stationary ergodic if the underlying probability space is actually an 
ergodic dynamical system with an action T and fix, uS) is such that fix +y,u>) = f(x, T y oj) for all x, y. 
In either case we may write f{x) instead of f(x,u>) and think of f(x) as a random function of x. 

Naturally, if / is now taking values in a finite dimensional space we will say it is stationary ergodic 
if all of its components in a given basis are stationary ergodic (in particular we may talk of stationary 
ergodic vector fields, matrix fields, and so on). 

Remark 2.1. Note that a function / is stationary ergodic if and only if there is some measurable 
function / : ft — > K (or — > E rf if / is vector valued) such that f(x,ui) = f(T x uj). 

Definition 2.3. A stationary ergodic domain is a map to — > 0{oj) C K d , to G ft such that 

(x,w) -> I ( w )(x) 

is a stationary ergodic function. Equivalently, 0{uj) is stationary ergodic if for some O G T we have 

O(w) = {a; | T^cj G O} 

The fundamental fact about ergodicity that we will repeatedly use is the ergodic theorem, see for 
instance [1! for a proof (in far greater generality that we use here). 

Theorem 2.2 (Ergodic Theorem). Given an ergodic dynamical system P-a.s. we have for any / G 
L" (ft), p G [1,00), that if f e (x,uj) := f(T x/e uj) then f e (x) E[/] in Lf oc {R d ) as e -> 0+. 

A particular and important case is when / is the characteristic function of a set O. 

Definition 2.4. Given a stationary ergodic domain O we have for any p G [1, oo) 

Ic.(aO :=P(0) inLf oc (M d ) P-a.s. 

We now introduce the function spaces we will be working with. 

Definition 2.5. For a given domain D we will denote by H{D) (iJo(-D)) the space of functions / G 
L 2 {D x ft) such that /(-,w) G i^O ) P-a.s. (f (-,<*}) G P-a.s.). Further, if a stationary ergodic 

domain is given then for any e > we introduce the set (which D and O are being used will be clear 
from context) 

B e := {(z,u) G D x ft | x G O e (cj) := e0(cj)} 

The set B e inherits the c-algebra and the product measure from D x ft, in particular we can talk about 
the spaces L p (B e ). Moreover, we define H(B e ) as the subset of L 2 (B e ) formed by functions / such 
that P-a.s. we have f{-,u) G H X {D n O e ), H (B e ) is defined as the closure of the set of / G H(B e ) 
which P-a.s. vanish in a neighborhood of 3D (but note they may be non-zero on dO e ). 

Definition 2.6. Given lu G ft and A C O e , we define the boundary of A with respect to 0(u>) as 

d u A := 8A\dO. 

As they are defined in different sets, we must define the convergence of functions / £ G L p (B e ). 
Definition 2.7. A sequence / £ G L P (B £ ) is said to converge to / G L P (D x ft), p G [1, +oo] if 

lrai ||/ e - /|Up ( b 6 ) = 0. 

Moreover, if p < oo we say / £ converges weakly to / (also denoted f e — /) if in the usual sense we have 

IbJ 6 / in L P (D x ft). 

It is also clear what we mean by convergence of functions G L P (B £( .) for a subsequence — > as well 
as what we mean when we say we have convergence in the P-a.s. sense. 



Quasistatic Droplet Percolation 7 

Definition 2.8. Given D C R d , e > and / e L 2 (B e ) we will say that v <E iJ(B c ) is a weak solution of 

^ ' 2.1 

= on dO„ 

if given any <f) £ H (Ti e ) we have 

/ Vw • V0 dxdP(u;) = f4> dxdP. 
Jb, Jb, 




Remark 2.3. Often, we will deal with (2.1) when D is a ball, thus it will be convenient to write 

B e r {x) := B r (x)r\O e . 
Definition 2.9. We will say v solves the Dirichlet problem 

-Aw =/ inDnO e , 
d n v = ondO e , (2.2) 
v = on a^D, 



if in addition to solving (2.1 1 in the weak sense we also have v & H (B e ). 



The mixed boundary problem (2.2 1 has a unique solution, which can be seen by applying the Lax- 
Milgram Theorem to a properly picked bilinear form in the Hilbert space i/ (B e ). 

2.2. Viscosity solutions of the Free Boundary Problem. Given a (deterministic) domain O C R d 
and a function u : O x M.+ — > R + , we introduce the positive phase or "droplet base" of u defined by 

D(u) := { u > } c0xl + . 

Let us also define D t (u) := {u(-,i) > 0} C O. On the other hand, any bounded set D C O determines 
a function uq € H 1 (C) by the conditions 

— Ait = 1 in Dq, d n u — on SO, u — in O \ D . 

Such Mo gives the droplet profile associated to the base Do- Then, given such a Do we consider the 
following Cauchy problem: to find a function u : O x R — s- M + with it(-, 0) = uo which solves 

-Am =1 in {u > 0}. 

<9 n u =0 on dO, (2.3) 

M t = |Vu| 2 on d u {u > 0}. 



It is well known that even for a smooth domain O and smooth initial data, solutions to ( 2.3 1 may become 
singular at some positive time, typically as the positive phase {u > 0} undergoes a topological change 
which induces a discontinuity in u. In this situation, some of the derivatives appearing (2.3) might cease 
to exist, and one must resort to a weak notion of solution in order to continue. Hence we now review 
the notion of viscosity solutions for (2.3), as it was developed in [2"9"j . 

We say xo £ dO is a regular point if dO is diffcrentiable at xq. Note that O may not have a boundary 
(for example O = M. d ). Also, E c O x R + will be called a cylindrical domain if it is of the form A x (a, b) 
for some open set A m O. 



Definition 2.10. A function tp G C{E) (E a cylindrical domain) is a classical subsolution of (2.3) in E 
if 



(i) if e CZ t {W >0}n£) and d{ip > 0} is locally in E. 

(ii) -A(p(x, t) < 1 in {tp > 0} H £. 

(iii) <y9„ < on regular points of dO n E. 

(iv) <y9 t < |V^| 2 on (d"{ip > 0}) n E. 
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Here Wp(x,t) is taken to be the limit of V<p(y,s) as (y, s) — » (x,t) from within the set {ip > 0}, and 
as before n denotes the outward normal of O. We say that ip is a strict classical subsolution if the 
inequalities in (iii) and (iv) are strict, lastly, classical supersolutions are defined by reversing inequalities 
(ii),(iii) and (iv). 

Definition 2.11. An upper (lower) semicontinuous function u in O X R+ will be called a viscosity 



subsolution (supersolution) of Problem (2.3) with initial data uq if u < uo {u > uq) for t = and if 



anytime we have a strict classical supersolution (subsolution) <fi of (2.3 1 in some cylindrical domain E 
with u < p (u > ip) on dpE then also u < tp (u > ip) in E. 

Given a locally bounded function u, we dehne its lower and upper semi-continuous envelopes: 

u(x,t) — liminf u(y,s), 

(y,a)->(x,t) ,^ 4) 

u(x,t) = limsup u(y,s). 

(y,s)-+(x,t) 



Definition 2.12. A locally bounded function iiinOx R + is a viscosity solution of (2.3) with initial 
data Uq if u is a viscosity subsolution and u is a viscosity supersolution, and u(., 0) = Uq. 



Now we can consider (2.3) for a (rescaled) stationary ergodic domain O e . 

Definition 2.13. Let us consider a stationary ergodic domain 0(lu), a compact subset Do C R d and e > 
0. A measurable function u e : B e x R + — > R will be called a viscosity subsolution(supersolution, solution) 



of ( f .1 ) with initial data Uq if P-almost surely the function u e (., w) is a viscosity subsolution(supersolution, 



solution) of (2.3) with O = O e (uj) and the corresponding initial data. 



For the case when O has a C 3 boundary we have a comparison principle, whose proof is parallel to the 
one in |29j . The regularity restriction on the boundary is to locally parametrize the Neumann boundary 
and use reflection argument to transform the flow in contact with the Neumann boundary into the flow 
without one (see section 3 of |30] where such paramctrization is performed in a similar context). We 
suspect the comparison to hold with less regularity, but we will not pursue this issue here. 

Definition 2.14. For nonnegative functions u,v : O — > K, which have compact supports and are 
continuous in the interior of their support, we write u -< v if u < v in {u > 0} and {u > 0} C {v > 0}. 
Following [29], we say u and v are strictly separated. 



Theorem 2.4. Let u and v be respectively viscosity sub- and supersolution of (2.3) with the domain 
O having C 3 boundary. Suppose u(-,0) -<! u(-,0) in O. Then u(-,t) -< v(-,t) for all t e [0,T). 

Let v be as given in the above theorem. When the initial free boundary d{v(x, 0) > 0} is C 1 , then 
Theorem 4.1 of [12] applies to v to yield that u(-,0) ~< v(-,e) for any e > 0. Using this observation and 
Theorem |2.4| we deduce the following: 

Theorem 2.5. Let u and v be as given above, and let d{v(x,0) > 0} be C 1 . If u(-,0) < v(-,0) then 
u(;t) < V*(;t) for all t £ [0,T). 



We do not clarify the uniqueness of viscosity solutions for (1.1) when dO e is less regular than C 3 or 



when d{v(x,0) > 0} is not C . In this case our result addresses the specific viscosity solution u £ given 



by the time derivative of the solution p e of an auxiliary obstacle problem (see (2.7) below) 



2.3. Main results. The main results require two assumptions on the sets O e , the second of which 
requires one further definition. 

Definition 2.15. An open set P C K d is said to be minimally smooth with constants (S,N,M) if wc 
may cover dP by a countable sequence of open sets {Ui}i such that 

(1) Each x € M. d is contained in at most TV of the open sets Ui. 

(2) For any x £ dP, the ball B$(x) is contained in at least one Ui. 
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(3) For any i, the portion of the boundary dP inside Ui agrees (in some cartesian system of coordi- 
nates) with the graph of a Lipschitz function whose Lipschitz semi-norm is at most M. 

We are finally ready to describe the two assumptions on O e . 

Assumption 1. The domains O e (cj) are given by 

O e {u)) := {ex | x G 0{u))} = {x \ T*u e O} 



where 0(uS) is a stationary ergodic domain and T the associated action (see Definition 2.3). 

Assumption 2. There are constants do,S,N, and M (independent of uj) such that P-a.s. the com- 
plement of the set 0(uS) consists of a countable union of bounded sets Pk(uj) (k £ N) such that we 
have 

d(Pk(u>), Pj(u>)) > do whenever k ^ j. 
Here each set Pfc(cj) is minimally smooth with constants (5,N,M) and has diameter smaller than d^ 1 . 

In light of Assumption [2j we make the following definition. 

Definition 2.16. A constant will be said to be universal if it is determined by the constants do, 5, N, M 
in Assumption [2] and the space dimension d>2. 

In a few words, Assumption [2] says that 0(u>) is almost surely obtained by populating WL d with many 
obstructions which are uniformly separated and sufficiently smooth (Lipschitz). In the homogenization 
literature and in the periodic case, this is known as a perforated domain. 

Initial Conditions. We take an initial droplet base which is a bounded open set Do C M. d and set 
Dq := Do H O e , we will assume that this set has no cracks, i.e. <9"Dq = <9 W D . Furthermore, we suppose 
that either Do is star shaped or that 8Dq is C 1 . This condition guarantees the uniqueness of the solution 
for the limiting problem (see [33]). 

With this notation in hand, the first of our two main results is the following. 



Theorem 2.6. Let O e satisfy Assumptions 1 and 2 and Dq be as above. Let u e solve (1.1 1 with initial 
positive phase Dq, and let u solve ( 1.2 ) with initial positive phase D (both in the viscosity sense), where 
/1 is as in Definition |2.4| and A as in Definition |7.2| Then the following is true: 

(1) The free boundaries d UJ {u e > 0} converge uniformly to d{u > 0} in the Hausdorff distance for 
t > 0. 

(2) Locally uniformly in t we have u e — > u P-a.s. in L p V p € [l,oo) (sec Definition 2.7). 

(3) Moreover, we have the following pointwisc limiting behavior 

u < lim inf u e < lim sup u e < u, 



where u and u are as defined in (2.4). 
(4) In particular, if u is continuous, then u e — > u P-a.s. in L°° locally uniformly with respect to t. 

As explained earlier in the section, even for smooth initial data the solutions u and u e may develop 
singularities, so classical solutions may not exist for all positive times and one must deal with viscosity 
solutions. We mention that topological changes for the positive phase of u are ruled out for its support, 
for example, when Do is star shaped (see Lemma 6.7). In any case, our result says that the oscillating 
free boundaries converge uniformly as e — > even if topological changes and discontinuities take place. 



The second result, which is needed for the proof of Theorem 2.6 deals with the uniform continuity 
of solutions to the linear elliptic problem (2.2). Since this is an interior regularity result we only need 



to state it for D = B r (x), and use the notation introduced in Remark 2.3 



Theorem 2.7. Suppose that O satisfies Assumption [2j There are universal constants C > and 
a 6 (0, 1) such that if v solves (2.2) in some ball B%, then 

IMIc(B; /a ) < C {r- a \\v\\ L * m + r 2 \\f\\ L ~) P-a.s. (2.5) 
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Moreover, if w > in £?^, we have a Harnack inequality, 



supv<C(inf v + r 2 ||/|| L oo) P-a.s. (2.6) 



2.4. Examples. Let us describe some of the random structures 0(ui) covered by Theorem 2.6 namely 
those for which Assumptions [l] and [2] hold (these are explained in Section 2.4). We also discuss cases 
that fall outside the scope of our current method as Assumption [2] fails (but "almost" holds) for them. 

Example 2.8. Site percolation with isolated obstacles: Consider a family of independent Bernoulli trials 
associated to each z £ Z d , defining a random subset G(w) C Z d where each z is included probability 
1 — p and excluded with probability p, independently of the others. Let Q = [0, 1]™ and I CC Q be a 
domain with a Lipschitz boundary, then define 

o(u) = m. d \ (J (x+z). 

The random set G(w) has a stationary distribution since the Bernoulli trials are identically distributed 
and independent, in fact G is ergodic so Assumption (JlJ is satisfied. 

Example 2.9. Other regular lattices A seemingly different example would be that of replacing Z 2 in the 
above example (when n = 2) with the triangular or hexagonal lattice, namely LcM 2 the set generated 
by linear combinations over Z of the vectors e l = (1,0) and e 2 = (1/2, 1/V2). We can run independent 
and identically distributed Bernoulli trials at each zgLto get a random set G(w) C L. Take a Lipschitz 
domain X that is compactly contained in the triangle defined by 0, e± and e2, then as before let 

o{Lj) = m d \ U (x+z). 

One could also consider a more complicated model where the domain I is itself random, as long as it 
is minimally smooth P-a.s. with respect to some constants (S, N, M) and it stays inside a fixed compact 
set inside the cell defined by 0, e\ and e%. Then O may look as given in Figure [l] (O c is given by the 
black region) 



X^X*X X*X ,-<^ x , X 1 
X'"V X^x^x' XX<#x""Xi 

; ¥ # X^X )XrX > 

x X 1 

Figure 1. A random domain O arising from percolation on a lattice. 



Example 2.10. (Modified) Irregular chessboard. Another example would be to consider an irregular 
chessboard with bounded cell sizes (compare with [151 Example 3.4]), except now instead of coloring a 
cell completely, we only color a rescaled cell half its size. The domain O obtained in this manner does 
not have a lattice structure, but it does however comply with Assumptions [T] and [2] 



Quasistatic Droplet Percolation 



11 



Example 2.11. General Site Percolation. If in Example |2.8| we do not assume that I CC Q then 
Assumption [2] no longer holds and our method fails. This is because we have no means of proving a 
global Sobolev embedding and Poincare inequalities if the perforations are allowed to be arbitrarily close 
to each other and are allowed to form large connected clusters. In fact, one does not expect the Sobolev 
and Poincare inequalities to hold P-a.s. in the general case. The best one can expect is that they hold 
in a given ball with a probability that approaches 1 very fast as e — > 0. This phenomenon is by now well 
understood in the context of discrete percolation, as discussed in Section 1.2. 

Example 2.12. Poisson Rain. For similar reasons to Example 2.11 we cannot handle the case where 
0(ll>) is built as the union of balls B ri (xi) where the points Xi are given by a Poisson point process and 
the radii r, are independent and identically distributed. Again, it would be worthwhile to investigate 
the extent to which the current techniques can be pushed to handle this case. This model is also known 
as the Boolean model in the stochastic geometry literature [231 Chapter 12]. 




Figure 2. An example where our assumptions do not hold is that of the Poisson rain. 



2.5. Strategy. In the analysis of ( 1.1 1 a constant challenge is the lower dimensional or singular nature 



of the free boundary condition. Namely, (1.1 ) can be written (in the sense of distributions) as 

dd{w>o} ~ Au e = I{„e >0 }, 
together with zero Neumann boundary conditions on dO e . The singular term is the time derivative of 
the characteristic function of the positivity set of u e . One obvious way to get rid of the singular time 
derivative is integrating the equation with respect to time. In particular, we may consider the function 

p e (x,t,ui) :— / u e (x,s,uj)ds 
Jo 

and seek an evolution equation for p e . This is particularly tractable in the one phase case, since u e is 
monotone in t. In fact, it was already known to Elliot- Janovsky [18] that whenever if is smooth enough 
then p e : O e x K + — > K is a non-negative function satisfying 

- Ap^x,t) = r(x,t)l {pC>Q} (x,t) 7 (2.7) 

and d n p e = on dO e , where 

r{x,t) = -l { ^ tY + ( I {r( , s )>o}(z) ds. (2.8) 
Jo 

This says that p e solves an obstacle problem |33j . The reason studying p e might be useful is that 
(assuming smoothness) it is not hard showing that 

{p e > 0} = K > o}. 
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Thus (heuristically) , one can study {u e > 0} using the methods available for the obstacle problem, which 
apply to {p e > 0}. 

In general, u e will not be smooth, and it may even become discontinuous at the onset of topological 
changes of {it e (-,t) > 0}, which may occur even if the initial data is smooth. This means we have to 



justify the above correspondence when u e is just a viscosity solution of (1.1). It is not clear how one 
may show that p e := J Q u e (-, s) ds is even a weak solution of (2.7) without differentiating u e directly. 

Following the method developed in |28j we invert the previous construction. Instead of defining p e as 
the integral of u e , we define it as a solution of the free boundary problem we expect it to solve (namely 
(2.7) and (2.8)). Only later we show that it agrees with the time integral of u e (see Theorem |6.1[). 



Definition 2.17. Let Dq C O e be as in Section 2.3 A measurable function p e : B £ x R + — » K will be 
called a weak solution of (2.7) if P-a.s. in uj we have that (x,t) — > p e (x,oj,t) belongs to the space 

A := {4> : O e (u) x R + -> R \ <p > 0, cj>{; t) G H l {O e {w)) V t}. (2.9) 

Moreover, with f e given by (|2.8[) we have 



Vp e - (p e - <j))f e dx > 



for any fixed time t and any <fi G A. 



The one to one correspondence of (1.1) to an obstacle problem suggests that to identify the limit of 
it e we should first homogenize the problem solved by p e . To do this, fix t > and let p = f Q u(-, s) ds, 



u being the solution to the homogenized problem (1.2) with initial data Do- The linear homogenization 
theory (which we revisit in Section [7]) says that if we can show that 

IoJ £ I {r >o} -» A^I{ P >o} P-a.s. in L 2 (R d ), 

to a p' defined in all of space and time, solving 



-a.s. 



(fi as in Definition 2.4) then p e converges 

- div(AVp') = /i ( -ji D 



Hp>o} 



Hp>o}- 



(2.10) 



Then, p and p' would solve the same elliptic equation with the same decay at infinity, and we would 
conclude that p' = p and therefore that p e — > p for each fixed time. 

Instead of proving ( |2.10[ ) directly, we will use the uniform continuity bounds for p e given by the theory 
in Section [3] and Section [4] These estimates give us enough compactness to extract a locally uniformly 
converging subsequence (in space and time) from any sequence p ek with efc — > (see Proposition 5.3). 



Then, all we have to show is that if p is the limit of such a subsequence, then it solves (for each fixed t) 

-div(AVp) 



/' 







(-Id§ ■+ 


%>o} ds^j 



Hp>o}- 



Remark 2.13. It is worth emphasizing how this last equation is different from (2.10). In the previous 
equation we had a function p' on the left and the positivity set of p on the other side, while here we have 
the same function p on both sides of the equation. 

In this case, since the problem has a unique solution for a fixed initial set Do, we would conclude 
that regardless of the initial sequence we have p — p, obtaining the homogenization for the obstacle 
problem ( |2.7| . The key to obtaining the above identity is the stability of the free boundaries d UJ {p e > 0} 
as e — > (Lemma 5.4), the proof of which relies on a uniform non-degenerate growth estimate for p e 
(Theorem |4.3[ ) . It is mostly for this single convergence result that we develop the elliptic theory in 
Section [3] and construct the special barrier in Section ffl 

Later in Section [6] we will show that any p e solving ( |l.l[ ) gives rise to a function u e which solves ( |1.1[ ) 
in the viscosity sense (see Theorem 6.1 1, and such that p e = J Q u e (-, s) ds. Using this Theorem and the 
strong convergence obtained for {p e > 0} we are able to prove that u e converges to the solution of the 
homogenized problem, proving Theorem |2.6| 
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3. Elliptic Regularity in O e 

3.1. Sobolev and Poincare inequalities. The goal of this section is showing that P-a.s. functions in 
O e satisfy inequalities similar to the standard Sobolev and Poincare inequalities in M. d . This is where 
the strength of Assumption [2] is most needed. 

We recall a classical result due to Calderon and Stein, concerning the extension of Sobolev functions 
defined in Lipschitz domains. For a proof see [JIJ Chapter VI Theorem 5]. 

Theorem 3.1. Let A C K d be a minimally smooth domain with constants (6,N,M), then there exists 
a bounded linear operator E : H 1 (A) —> with its operator norm depending only on 6,n,M. 

This extension theorem will allow us to exploit the desired inequalities from M. d and pass them to O e . 

Proposition 3.2. Let A C R d be a connected domain which is minimally smooth with constants 
(S,N,M) and diameter no larger than (Lq . There is a constant C = C(S,N,M,do) such that for any 
6 G H 1 (A) we have the following Poincare type inequality 



10- {4>)a\\l*(A) < C||V<?!>| I £2 (A) 

where (4>)a '■= / </> dx. 

J A 



Proof. We argue by contradiction, relying on compactness. Suppose that for some combination of 
(5, N, M) and do such a constant C does not exist. Then for each j € N we can find a connected, 
minimally smooth domain Aj with respect to (6, N, M), with diameter no larger than d$ , and a function 
cj>j € H 1 (Aj) such that 

Hi ~ (^)a,||l*(A,) > V\\V^\\ L 2 (Aj) > 0. 
Moreover, after a translation, we may assume that all the Aj's contain the origin, so that Aj C B 2d -i (0). 
On the other hand, each <pj is not identically equal to {<f)Aj in Aj, so we may introduce normalized 
functions 

$3 ■= Uj - (^AjWjt^ifo ~ i^Aj), 

so that 

(^•)^=0, llV'jlU 1 ^) = 1) and II V-0il| J L2 (Aj . ) < 2~K 
All the sets Aj are minimally smooth with the same constants, then Theorem 3.1 says we can find 
constant Co = Cq(S, N, M) and functions -0j G with jpj = ipj in Aj, such that 

\$j\\ H im < Co, Uj\\m[Aj) < c (i + 2- j ). 

Since ipj extends ipj, we have HV'jIIl 2 ^) ^ 1 f° r each j. Moreover, given the gradient estimate we may 
assume, by passing to a subsequence, that ipj converges strongly in L 2 (M. d ) and weakly in H 1 (M. d ) (and 
then in fact, strongly in H 1 (R d )) to some Voo € H 1 ^). 

On the other hand, the fact that Aj C B 2d -i for any j and that is minimally smooth show that after 
passing to a subsequence Aj converges in the Hausdorff topology to some set A^. Furthermore, it is 
not hard to see that A^ will be connected and minimally smooth with respect to (8,N,M), and in 
particular it has a non-empty interior. In conclusion, for every j we have 

$j) Aj = o, II^-IU^.) = i, ||v^|U 3(Aj) < 2-j, 

from where it follows that (V'oo)a 00 = and Vi/'oo = in Aoo, and since Aoo is connected ipoo must be 
identically zero a.e. in A^, but we also must have ||V'oo||l 2 (.a 00 ) = 1 which is impossible. This contradicts 
the existence of the original sequence <j)j and proves the proposition. □ 

Although Poincare's inequality itself is used in proving elliptic estimates, here we need it only to 
prove the Sobolev inequality, which is what will we use in the proof of Theorem |2.7| We also recall the 



definition of a universal constant given in Definition 2.16 which will be used extensively. 
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Proposition 3.3. Assume O satisfies Assumption [2j then there is a universal C such that for d > 3 

\\<t>\\ L *'(Oc)<c\\v<f>\\ L z m , v^if'taH), (3.i) 

for P-almost every u>, here 2* := 2d/ (d — 2). 

Remark 3.4. As it is well known in the case d = 2 and O = M 2 the above estimate with L°° on the 
right hand side does not hold. Instead, the relevant estimate for d = 2 says that u is a function of 
bounded mean oscillation. Namely, there is a universal constant C such that for any x and r > 



|u - (")x,r| dx < C\\Vu\\ L 2(o,)r , (u) x ,r ■= / u(y) dy. 

Although all of our estimates hold for d = 2, as long as one uses the above inequality in place of the 
Sobolev embedding. For the sake of presentation we will only work our arguments explicitly when d > 3. 

Proof. For any k,ui, we denote by P k (u>) a do/4-neighborhood of P k (uj). By assumption, almost surely 
the sets P k (u>) are minimally smooth (with constants independent of k,uS). Then, by Theorem 3.1 case 
we know that there is an extension operator E k such that 

E k : H\P k (u) \ PfeH) -> H\P k (Lj)) 

where for some C = C(S, N, M) (and thus independent of k), we have 

\\ E k<l>\\imp h (u)) - C W^WliHP k (u)\P k (u)y 
We define new extensions E k : i/ 1 (P fc (^) \ -Pfc(^)) -t H 1 {P k {w)) by 

E k( j, := - (<fi) k ) + 
and, putting them all together, an extension E : i? 1 (C) — > H 1 (M. d ) given by: 

!4>(x) whenever x G O 
E k (f>(x) whenever x £ P k (u)). 

Now, in P k (oj) \ P k (w) we have E k <j) = ((f)- (4>) k ) + ((f)) k = 4>, thus 

\\VE k <f)\\ L 2 { p k(uj)) = \\VE k ((f>- (</>)k)\\L*(,p k (w)) < C\\<f>- (<f>)k\\ H i(P h ( u )y 



E(f)(x) 



Then, due to Assumption [2] we may apply Proposition 3.2 to P k (ui), so there is a universal C such that 



\\<t>- (4>)k\\ m (p k { u )) < C||V^|| £2(j p fc(w)) 
since this holds for every fe, we have proved that 



[ \VE(f)\ 2 dx<C 2 [ \Vcf)\ 2 dx. 



Therefore, 

f \\7Ecf)\ 2 dx < (1 + C 2 ) f I V0| 2 dx. 

JWL d JO 

Next, observe that, since E(f> = <f> in O we get the bound 

IML»-(O0 < ll^ll^(R^)- (3-2) 
Applying the standard Sobolev embedding we get 

II^IL 2 *(Oi) < Cn||VS</>||L2( R d), 



so, using the bound (3.2 1, we obtain for a universal C 

This is the desired estimate when e = 1. These norms scale in the same way, so we may use the change 
of variables x — > ex to obtain the estimate for general e > by reducing it to the inequality above. □ 
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Next, we prove an inequality similar to Poincare's inequality, except now we assume the function 
vanishes in a fixed portion of the ball. It extends a corresponding inequality in R d used by Moser in his 
original proof of the parabolic Harnack inequality |38) . It seems weaker than De Giorgi's isoperimetric 
inequality for H 1 functions, which says that if u G H 1 (Bi) and < u < 1, then 



|{u = 0}||{ U = 1}| 2 /" < C n \\Vu\\ L2iBl) \{0 <u<l}\ 



1/2 



This inequality was used by De Giorgi in his proof of the Holder regularity for elliptic equations |16| . 
We were not able to prove an analogue to De Giorgi's inequality due to the irregular geometry of the 
perforated domain. However, we observe that for the purposes of regularity and Harnack's inequality, 
our extension of Moser's estimate is enough. We recall the statement of Moser's estimate, see [3H Lemma 
2] for a proof. 

Lemma 3.5. Let (f> G H l (B r ) and d > 3. Then for any m G (0,1) there is a positive constant C 
determined by d and m such that 

f - f 



W dx ) < C k / |V0r dx+ \ 4> A dx 



where A is any subset of B r such that \A\ > m|S r |. 
The extension of Lemma l3~5l is as follows. 



Proposition 3.6. Let cf> G H 1 {B*{xq)) and d > 3. Then for any me (0, 1) and d > 3 there are positive 
constants C and c\ determined by 5, N, M,do,d and m, such that P-a.s. the lower estimate 

\{x G B e r (x ) | 4> = 0}| > m|B^(^o)| 

guarantees the bound 

\\<t>\\L2(B' cir (x )) < Cr 2 1| V0|| i» (B»(xo))- 

Proof. As in the previous proposition, let us consider first the case e = 1. Let I C N be the set of indices 

Consider the domain 



k such that P k (uj) n ^ 0, and let Pkioj) be as in the proof of Proposition 



3.3 



fee/ 



Arguing as before, there exists a bounded extension operator i?o : -ff 1 (J r ) — > whose norm depends 

only on O(oj). First suppose that r > 1 and observe that in this case J 7 C -B<> for some universal C . 
Since we have a function E <f> which vanishes in a portion of Ben Lemma |3.5| says that 



U^Ul^) < C'r 2 \\VE Q cp\\ LHBcr) , 
but now one can see that for larger (but still universal) constants C and C" , we have 

||v^^|U a(Bor) < cViiv^ii^^). 

Hence it follows that as long as r > 1 we have 

UWlhbv < C"r 2 \\V<f>\\ L 2 (Bhiir) . 

If we can prove the inequality for r < cq for some universal Co then we would have proved the proposition 
via a doubling argument. Let cq be small enough so that r < cq guarantees that B r (xo) (1 d£l is the 
graph of a Lipschitz function in some Cartesian system of coordinates. In this case, a Lipschitz change 
of variables allows us to reduce the estimate to the case of a half ball. Here we may argue by a reflection 
argument to reduce it again to Lemma [33J and the proposition is proved. 

□ 
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The proof of Theorem |2.7| has two stages, first for / = we prove Harnack's inequality using a 
combination of De Giorgi's and Moser's approaches (to get the L°° estimate and oscillation estimate, 
respectively) . 

Once this is done, we prove a Stampacchia-type estimate to handle / ^ 0, with / G LP (p > d/2). 
Combining this estimate with Harnack's inequality from the homogeneous case we get the general result. 
For a thorough presentation of similar ideas for the classical case of a uniformly elliptic operator with 
bounded measurable coefficients in K d , see [33J Chapter II]. 

3.2. Harnack inequality (homogeneous problem). As we start the proof of the regularity estimates, 
we recall the notion of a subsolution. 



Definition 3.1. A function v is said to be a weak subsolution of (2.2) if for any positive test function 
cj) whose support is compactly contained in B r , we have 



Vw • V4> dx < / /</> dx 

Wc define a weak supcrsolution similarly by reversing the last inequality above. Note that a function 
which is both a subsolution and supersolution in this sense is a weak solution (see Definition 2.8). 

Remark 3.7. Given w € LL 1 , the function w\ := (w — A)+ is also in LL 1 . It is well known and not hard 
to see that whenever w is a weak solution to (2.1 ) then w\ is a weak subsolution for any Ael. 

We next recall the Cacciopoli or energy inequality for subsolutions. 

Proposition 3.8. Let w be a weak subsolution with / = 0. Then for any A € K and any Lipschitz 
function rj with compact support in B r we have the inequality 

/ \V( V w x )\ 2 dx< I \V V \ 2 (w x ) 2 dx (3.3) 
Jo c Jo, 

Proof. The function w is a weak subsolution, so for any positive 4> & Hq(B t ) we have 

Viu ■ dx < 

o € 

Putting cf> ~ r/ 2 w\, we get 

/ V«i • (2r]wWr) + ^ 2 Vw A ) dx < 
we "complete the square" on the left by adding rj 2 \Vri\ 2 w\ to both sides, getting 



Vw • (277w A V?7 + ?? 2 Vu;a) + (wa) 2 |V?7| 2 dx < / \Vr]\ 2 (w x ) 2 dx. 
Indeed, since V(rjw\) = r)\7w\ + waVjj the integral on the left equals 

\V(r]Wx)\ 2 dx 

which gives the first inequality. □ 

Before we prove the first L°° estimate, we need the following elementary proposition. 
Proposition 3.9. Let be a sequence of positive real numbers and a, ft, S > numbers such that 

A k+1 < a2 l3k A 1 h +s for each k e N. 



Then lim A k = if 

k— >-\-oo 



An < 2-^V 1 / 5 . 



Quasistatic Droplet Percolation 17 

Proof. We shall show that A k < 2~ kfM Ao for some /i > 0, from where the proposition will follow. When 
k = this is trivially true, so suppose that for some k^ we have 

Then, by the hypothesis 

A ko+1 < a2 fSko A 1 k + a s < a 2^ k «2- ko ^ 1+ ^A 1 +s , 

so we are done as long as 

If we set fj, = P/ 5, this last inequality is equivalent to 

^<5 ^ 2^(' £ n + l)M-^ fe o+fcoA 1 (l+i5) Q ,-l _ 2^t J -^0 k o+^Sk o ^-1 _ 2 _ ' 1 q : ~ 1 

so Ao < 2~P/ & a" 1 /* 5 guarantees that Afe -> as we wanted. □ 

We are now ready to prove that subsolutions are P-a.s. bounded from above, at least when / = 0. 
The proof relies on the fact that if v is a subsolution then v\ satisfies both the Sobolev and energy 
inequalities for all Ael. 



Lemma 3.10. Let v be a weak subsolution of (2.2). Then for some universal C we have 

C f 

sup v 2 < — - / v 2 dx P-a.s. 



Proof. As it is usual in the energy method, for each k £ N we define 

r k :=r(l + 2- fe ), 
M k :=Mo(l-2~ fe ), 



and subsolutions (see Remark 3.7) v k := (i> — M k ) + . We also define cut-off functions r\ k £ Lip(_B r ) 
where < 77^ < 1 everywhere, rj k = in R d \ B rk _ 1 , r\ k = 1 in B rk , and |V7?fc| < 2 k jr. Now, let 

A k ■= / (rj k v k ) 2 dx 



we shall see how large Mo ought to be to guarantee that 

lim A k = 0, 

k— >-+oo 

in which case we may conclude that u < M almost everywhere in B r . 
Applying Holder's inequality, we get 

A k <( [ {r)kv k ) 2 *dx) 2 \{ Vk v k > Ojf-l 



Then thanks to the Sobolev type inequality (3.1) and the Energy inequality (3.3), 

A k <C J \V(r] k v k )\ 2 dx\{r] k v k >0}|5 < C J \Vrik\ 2 v 2 k dx\{r} k v k >0}|5 
since \Vr] k \ < 2 k /r, we get 

A k < C2 2k r- 2 f v 2 k dx\{ Vk v k >0}\% 

J {Vk^0} 

Also from the definition of r\ k we have that 77^—1 = 1 whenever r\ k ^ 0. Moreover, w k =/= means that 
w > Mfe, and since M k+ \ = M k + M 2^ < - k+1 \ these two observations imply that T] k -iv k -i > Af 2 _ ( fe+1 ' 
in \jj k v k 7^ 0}. This gives the measure estimate 



\{VkV k > 0}| < 2 2 ( fc+1 )A/ - 2 / (vk-iVk^) 2 dx 
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and we finally get 



From here, applying Proposition 
lim Ak > it is necessary that 



3.9 



4 

with a — C$r~ 2 (AMq 1 ) d } f3 = 6 and 5 — | we conclude that for 



A < 2~ 6/s2 (Ct- 2 A 2S Mq 2S ) s = 2- 6/s2 C~Sr6 4 2 Ml 

In other words, 

M 2 > 2- 6/s2 ch-'h- 2 A . 

This shows that 

2 C J/2 C f 9i 

sup v < —jAJ < — 7 / v dx. 

□ 

Likewise, if v is a supersolution, we conclude it is bounded below since — v is a subsolution. This shows 
that weak solutions of (2.2 1 are in L^ c . Next, we prove a very weak version of the Harnack inequality. 



Lemma 3.11. Suppose that v is a supersolution of (2.2 1 in B^, and that v > 0. If M > is such that 

\{x € B e r | v{x) > M}\ > m\B% 

then we have 

inf v > c M. 

where cq depends only on O and m and c\ is a small universal constant. 

Proof. We start by deriving an energy inequality similar to the one used in the upper bound. Fix 7 > 
(later we will let it go to zero), taking (f> = r/ 2 (v + 7)~ 1 as a test function in the supersolution equation 
for v, we obtain 

0< J^v [-0 + 7)~VVv + (t> + 7)~ 1 277V?7] dx. 
Moving the first term to the left, this becomes 

.,2/ 1 _,\— 2|V7,.|2 „J ^ n I „/„. 1 \— 1 



ff{v + i) \Vv\ t dx<2\ r)(v + j) L Vv-Vr]dx. 
Then, using the Cauchy-Schwartz inequality on the the right hand side and rearranging, we obtain 

?7 2 |(« + 7r 1 V(v + 7)| 2 dx < 4 / |V?7| 2 dx. 



We recognize on the left hand side the gradient of log(i; + 7), so, taking 77 to be a smooth function which 
is identically 1 in B r / 2 and vanishes outside B r we arrive at the estimate 



/ lVlog^ + 7)! 2 ^ <C d r- 2 \B% 

J K/2 

Next, consider the function v = max{— log(^p^), 0}, and note that 

/ \Vv\ 2 dx< / |Vlog(a; + 7)| 2 dx. 

Take 7 small enough so that v + 7 > M in a set of measure greater than 2 _1 m|S^|, then we can apply 
Proposition |3.6| and conclude that 

f v 2 dx<Cr 2 [ \\7v\ 2 dx < C\B% C = C{d,m). 
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On the other hand, it can be checked that v is a subsolution. Hence, Lemma [3. 1Q| yields that 

sup v 2 < ^C\Bl\<C. 



B« ir/4 r 



Where C is again a universal constant. Using the definition of v this says that 

-log[0 + 7)/M] <CC^ [v + i)/M > e~ c . 
Letting 7 — > 0, we conclude that 

inf v > e~°M, 

which proves the lemma with cq = e~ c , so that Co is a universal constant. □ 
Corollary 3.12. There are universal constants C, c\ > such that P-a.s. we have 

sup v < C inf v 



for any v which is a nonnegative solution of (2.2 1 in B%. 
Proof. Indeed, if v is a solution, it is bounded from above, so let 

M = \ sup?;. 



Applying Lemma 3.10 to (y — M) + , there exists a constant C such that 

MV < C f (v - M) 2 + dx < CM 2 \{x e B e r I v > M}\ 

and we conclude that 

\{xeB^.\v> M}\ > G- x \Bl\. 
This, together with the fact that w > 0, allows us to apply Lemma [3.11| which gives the bound 

c M < inf v. 

~ Bi ir 

Putting this together with the definition of M the desired inequality is obtained. □ 
Corollary 3.13. There are universal constants /i,ci € (0, 1) such that if v is a solution in B^ r , then 

osc f < /i osc ?; P-a.s. 

In particular, there are a = a(O) £ (0, 1) and C — C{0) > such that 

\HcHBi ir ) < M v W^m P-a.s. 
Proof. If v is constant in B f Cir there is nothing to prove, so let M > m be given by 

M = sup v, m — inf v. 

B>„ B< 



Then v* = j^Zm ^ s harmonic and positive in B r U O e . If c\ < 1/2 then Corollary 3.12 yields 

sup v* < sup v* < C inf v* . 

B^ r B* /2 B^ r 

Suppose first that C inf v* < |, then (note that oscu* = 1 by construction) 

osc v* < sup v* < h — s osc v* . 

B kr B kr 2 2 B% 

Otherwise, we have inf v* > hC~ 1 , in which case 

Ft" ~ 2 

"r/2 



OSC V 



< 1 - inf v* < 1- \C- Y = (1 - \C- x )omv* 
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In either case, with fj, = max{i, 1 — \C 1 } € (0, 1) we have 

osc v* < ji osc v* 



this inequality for v* immediately implies the desired inequality for w. □ 

3.3. Stampacchia's Maximum Principle. The last tool we need in order to prove Theorem |2.7| is 
Stampacchia's maximum principle, which deals with / ^ and zero Dirichlet boundary conditions. 



Lemma 3.14. We have P-a.s. that for any v € Hq solving (1.2 1 in B% and p > d/2 (d > 3), 

2-^ 

IMU»(B«) < Cor P\\f\\ L p(B' r )- 
Proof. The proof is a variation on the theme seen in Lemma |3.10| Fix M > and define for each fceN 



Vk 



M(l - 2- k ) 
(v - M k )+ 

\{v k > 0}| = \{v > M k }\. 



Since 



\{(v-M) + >0}\ 



lim A k 

k— >-\-oo 



we will conclude that v < M a.e. in B e r by showing A k goes to zero. As before, recall that anywhere 
where w k > we also have v k -i > 2~ k M, so by dividing and multiplying v k -i by 2~ k M we get the 
estimate 



A k < (2 k M- 



-1\2* 



vi_ 1 dx<2* k M-> C s ||V« fc _i| 



(3.4) 



Here we used the Sobolev embedding (3.1) for the last inequality. Now, since v k -\ vanishes outside B% 
we can use the equation and Holder's inequality to get 



|Vufc. 



dx 



\v k f \ dx < ||«fc-i||2«||/lk 2 < i /(< i + 2 )(A fc _ 1 )- 



l{v k >0} 

Using (3.1) again and dividing by ||Vufc_i||2 on both sides we have 

||VUfe_i|| 2 < Cs||/|| 2 d/(d+2)- 
Since the integral for g is over A k -i we can use Holder's inequality again to get an extra term, 

||V« k _i|| a < C s || 5 || 2 < Cs||fl|| P |{«*-i > 0}|^N = C s \\f\\ p A^ p . 



(3.5) 



Combining (3.4) with (3.5) we prove the relation 



A k <2 2 * k M- 2 *Cf\\g\\£A k _\ 



d+2 _ 1 
2d p 



Then, we can apply Proposition 3.9 to the sequence A k with 



* = 2*(^-i)-l = ^(4-f)>0 



to conclude that lim Ak — whenever M satisfies the inequality 

fc— >-+oo 



A a <Co\\f\\ 



-2*8~ 



As in Lemma (3.10) this shows v < M for some M for which the reverse inequality holds, namely 



M 2 * r 



< Co||/||f ^"Ao- 



Since 5/2* = ^ (4 - = s( 2 ~ p)> we conclude that 



M<C \\g\\ P A» 



Quasistatic Droplet Percolation 



21 



Now, \A \ < \Bi\r d , therefore A^ 2 < C n r dS / 2 * — Car 2 d and the lemma is proved. 
Finally, we put all the pieces together for the full elliptic estimate. 



□ 



Proof of Theorem 2.7 For any p < r/2 we seek to bound 
In order to do this, let us decompose v as 



osc v. 

B p {x ) 



v = Vq + V\ 



where vq = v on (dB p (xo)) D O t and solves (2.1) with / = 0. Then 

Applying Lemma |3.14| with p = +oo 

o 

Bp 

and by Corollary |3. 13| and Lemma |3.10| 



osc v < osc v + osc Vi . 

B p (x ) B p (x ) Bp(x ) 



osc wi < 2 sup \\vi\\ < Cp 2 \\g\\ 

B P {x ) Bp(x ) 



Therefore 



osc v < C\\v\\ L 2p a 

Bp(x ) 



osc v< C(M\l> + p 2 - a \\9\U) P° 

Bp(x ) 



and the first part of the theorem is proved. Now suppose that v > in B!p, then if vq and Vi are as 



before, it is standard to see in this case that vq > everywhere in B e r , so by Corollary 3.12 



sup fo < Co inf w • 

K,2 B '/2 



Therefore 

sup v < Cr 

Moreover, inf vq < inf v + \\vi\\l^^ and thus 



sup v < Cq inf vq + sup Vi 



■V/3 



r/2 



sup v < Cq inf v + 2||i>i|| i oo( B « 



r/2 



r/2 



and using Lemma 3.14 again we finish the proof. 



□ 



4. Weak free boundary results 



We are now ready to anal yze the behavior of the free boundary <9"{p e > 0} for the auxiliary problem 

Observe that if p e was uniformly C 1 in {p e > 0} and had a non- 



12. 7) introduced in Section 



2.5 



zero derivative on d ul {p' i > 0} then the implicit function theorem would say that 5 w {p e > 0} is a C 1 
hypersurface. The problem is that this is never the case for the obstacle problem, since p € vanishes 



quadratically at its free boundary (see Lemma 4.2) 



On the other hand, we shall show p e (-,t) vanishes not faster than quadratically (Theorem 4.3). More 



concretely, for each fixed t the supremum of p e (-,i) in a ball centered at a free boundary point is 
comparable to the square of the radius of that ball, and most importantly, this holds uniformly in e. 

Although this is far from a C 1 estimate on the free boundary, this non-degeneracy estimate is enough 
to prevent the free boundary from developing many thin fingers as e goes to zero. This will guarantee 
the stability of {p e > 0} as e — > 0, that is, if p e converges along some subsequence to a function p, then 
{p e > 0} converges to {p > 0} along that same subsequence (see Lemma 5.4). 
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4.1. Estimates with respect to x. First of all, we remark that solutions to the obstacle problem in 
O e in the sense of Definition |2.17| are continuous in x for each single t. This follows by extending a 
classical theorem in potential theory to O e . The proof does not require uniform estimates in e, but using 
directly the Harnack inequality and Stampacchia maximum principle from Section [3] one may adapt the 
proof from the Euclidean case. 

Lemma 4.1. (see [7J Theorem 1]) Let O be a deterministic domain satisfying Assumption [2] and p a 
nonnegative function in iJ^O) such that for some / 6 L°°(0) we have 

-Ap = f in {p > 0}, d n p = on {p > 0} n dO 

and — Ap > / in O, d n p > on dO. Then p is a continuous function in O. 



From this qualitative continuity lemma, we may use Theorem 2.7 (again) to show p e vanishes quadrat- 
ically at its free boundary, arguing as in the standard obstacle problem Theorem 2]. 

Lemma 4.2. Suppose p £ is a solution in the sense of Definition |2.17| Then, for any e, t > we have 

-Ap e = f e l {r>a} in O e x {*}, d nP e = on dO e x {t}. 

In particular, p e (-,t) is Holder continuous uniformly in e, and for any x £ d"{p e > 0} we have 

sup p £ < q|/ e (-,t)|Ur 2 . 

B$(x)x{t} 

where C is a universal constant. 

Proof. Fix t (and we will simply write p e instead of p e (-,t)), and let x G d u {p e > 0} and r > also be 
fixed. Recall that — Ap e = f e in the interior of {p e > 0} and that in general 

-Ap e > f in O e , d n p e > on dO t . 

Since we already have the validity of the equation away d ul {p t {-, t) > 0} we only need to focus on x near 
d UJ {p t (-,t) > 0}. In particular, we may assume without without loss of generality that f e (-,t) < in 
B e r {x), once we have obtained the equation in this case, the bound on supp e would follow for all other r. 
Then, we decompose p e as p e = p\ + p\, where 

-Ap\ = f e in B e r (x), d n pl = on dO e , p\ = p e on d^B^x), and 

-Apl > in B%{x), d n p\ > on dO e , p e 2 = on d u B e r {x). 

By the maximum principle and the fact that / £ < we see that p\ < p e in B%. Moreover, since p\ > 
in dB e T we may use Lemma 3.14 to obtain an opposite bound, 

pi > -cwnur 2 

for some universal constant C. In this case, the function v = p e + C||/ e || 00 r 2 is nonnegative and solves 

-Av = f e in B e r (x), d n v = on &*B e r (x). 
Therefore, applying Harnack's inequality (Theorem |2.7p we get 

sup v < C (v(x) + ||/lccr 2 ) . 

From Lemma 4.1 we know p e is continuous and that x £ d ul {p t > 0}, so we have p\(x) < p e {x) = 0. 
Combining this with v(x) = p\(x) + C||/ £ || 00 r 2 and the Harnack inequality above we obtain 

sup p\<C\\r\\r 2 

for another (still universal) constant C > 0. Since p\ > and it vanishes on d^Bf, the maximum 
principle says that the supremum of p| must be achieved in the closure of the set p e = 0, and there we 
have p\ = —p\ so that p\ < Cr 2 . 

In conclusion, for every x € d u> {p t > 0} and any r > we have 

< p e < Cr 2 in B £ r (x). 
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In particular, p e vanishes quadratically on d UJ {p £ > 0}, from where it follows (by a standard test function 
argument) that d n p e = on dO e n d u {p" > 0} and that -Ap e = I{ p «>o}/ e in O e . 
In this case, the Harnack inequality from Theorem |2 . 7| applied to p e > yields 



sup p e <C[ inf p e +r 2 \\f e (-,t)\\ 

and since p e (x,t) — 0, the infimum above must be zero, and the lemma follows. □ 
Now we focus on showing that p e (-,t) does not vanish faster than quadratically. 



Theorem 4.3. There is a universal Cq > such that if a; € {p e > 0} and r < |c?(x,Dg), then 

sup p e (-,i) > e~ Co(t+1) r 2 

The argument in [7], which is later used in [28] fails under the presence of the random domain O e , as 
it relies on the fact that the function ^|a;| 2 is a supersolution and thus can be used as a barrier. This 
is no longer the case for our current problem case since <^|a;| 2 does not always have the right Neumann 
boundary data on dO t . 

Therefore, we need to come up with a replacement, which one may think of as a corrected paraboloid 
for the domain O e . The construction of such barrier is the the content of the following key lemma. 

Lemma 4.4. There is a universal constant a such that P-a.s. for any e and any xq £ O t there is a 
continuous function v : O e — > K, with v £ Hl oc (O c ) and such that (in the H 1 sense) 

-Av > -1 in O e (uj) 
d n v > on dO e {uj) 

v > ar 2 in O e \ B^(xo) for any r > 
v — for X — Xq. 

Remark 4.5. The construction of v is not so straightforward and will take several preliminary lemmas. 
It is worth remarking that, one cannot simply invoke known homogenization results [■¥.)[ 142) and use 
strong L 2 convergence for elliptic equations in perforated domains and borrow a barrier from the limiting 
equation. This approach will not suffice, since later arguments use strongly both that v — at xq and 
that it grows quadratically away from it, and such a pointwise condition is not stable under perturbations. 

The first thing we will need is a basic geometric property of O e , which follows from Assumption [2j 

Proposition 4.6. There are constants L > I > such that P-a.s. we have 

there exists y z £ O x s.t. Bi{y z ) C O x n Ql(z) V z £ Z d , 

where Ql(z) denotes the cube {x £ R d | \x — Lz]^ < L}. 

Proof. Let Pfe(w) denote the perforations giving rise to C(w) and do,8,N,M the constants related to 
the regularity of O, as described in Assumption [2] (Section [2| , also let L be the smallest positive integer 
larger than 2d^ 1 . For any cube Q, we will denote by Q* the cube with same center and half the diameter. 

Let us find for each z £ 1 d Then, for each z £ % d the cube Ql(z) must be such that Ql( z )* lies 
entirely in O(oj). In that case we may take y z = Lz so that BL(y z ) c 0(lu). If not, there is some point 
x and some k such that 

x£(dP k (w))nQ L {z) 

In particular, since Pfe(w) is a minimally smooth domain, there will be a ball of radius r$ = ro(S, N, M, do) 
in some neighborhood of x which will lie entirely inside 0(lj), then in this case take y z as the center of 
this ball, and let I = min(r, L), and the proposition is proved. □ 



We shall now use the family of balls {Bi(y z )} zeZ d from Proposition 4.6 and a standard Whitney 
extension [JT] to construct the barrier far away from xq (see Figure [4]). 
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Figure 3. The geometric configuration considered in Proposition 4.6 



Lemma 4.7. There are universal constants do, fli, with ao £ (0, 1) ao < a\ and ro > such that P-a.s. 
for any e and any xq £ O e there is a continuous function v : O e — > R, u € i// oc (C £ ) and such that 



|A«| 
<9„w 

U 
U 



< 1 

= 

> aor 2 

< air 2 



in O £ (,)\B; £ ( l0 ), 
on dO e (u), 

in O e \ 23' (xo) for any r > roe, 
in O e \ B e r {xo) for any r > r e. 



The first two equations are understood in the H 1 sense. 

Proof. We start by noting that it is enough to build such a barrier for e = 1 and an arbitrary point 
e Xo £ 0\. For if v satisfies all of the above properties for e = 1 and the point e~ 1 Xo, then the function 
e 2 v(x/e) satisfies all the properties listed in the statement of the lemma. 

Let I, L, {y z } z and Q z :— Ql{z) be as in Proposition 
all of K d , and we have the inclusions 



4.6 



The non-overlapping cubes {Q z } Z £Z d cover 



%)cOinft, Vzez d . 

Given x £ M. d we let y Zx denote the closest point to x among all the {y z } z - Take R > and define 
Ir := {z £ Z d | Bi(y z ) c B R (xo), \Lz — aso | > !}• Then, we introduce the domain 

D« = Bji(so)\ (J B l/4 (y z ). 

zeln 

Note that Di? is non-empty for R large enough. We will build approximate functions vr in Dr with 
controlled growth in R to extract a (locally) converging subsequence as R — > +oo. For any large enough 
R we define vr in as the solution to the following obstacle problem 



ur 

O n VR 

|Vur| 
vr 



> in Dr, 

= 1 in {vr > 0}, 

= ondO 1 UdB R (x ), 

= on > 0}, 

x \ 2 on 9 J B ;/4 (y z ), z £ Ir. 



\Lz- 
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This is understood in a weak sense modifying Definition 2.17 accordingly (the existence of vr is done 
following Appendix |A| . Then, arguing as in Lemma 4.2 one can show that vr satisfies 



Av R = I{ VR>0 } 

9 n VR = 



in D R , 

on <90iU dB R (x ). 



On the other hand, using a barrier argument, one can show that 

v R > in B C() i(y z ), V z £ Ir 

for some small (universal) cq > 0. In this case, when z £ Ir, since vr solves an elliptic equation in 
a smooth domain with smooth coefficients we have that vr is in C 2 (B Co i/ 2 (y z ) \ B C[) i/4(y z )) and with 
universally bounded second derivatives (that is, bounded independently of R) Thus, invoking again an 




Figure 4. A schematic description of the barrier v defined in O. 



extension operator [41] we extend vr it to a C 2 function in all of Bi/ 2 (y z ) (for each z £ 7r), resulting 
in a function defined in all of B r (xq), resembling the one in the figure above. Let us keep calling this 
extension vr. Then, for some universal C\ > we have 

\Av R \ <Ci mB R (x ), 
d n VR = on dO i U dB R (x ), 

vr = \Lz — xo\ 2 somewhere in Biu(y z ), V z £ Ir. 

Next, note that r\ := \pnL is such that if x £ Dr and B\ x (x) C Dr then we have 

sup v R > \Lz x - x \ 2 , 

(x) 

where z x £ Z d is defined as that for which \Lz — x\ is the smallest. From here, using the Harnack 
inequality from Theorem |2.7| and the definition of z x , we can see that 



inf v R >C-\ sup VR-r 2 ) > C-\\Lz x - x \ 2 - nL 2 ). 
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Moreover, since z x was picked so that \Lz x — x| L we have (by the triangle inequality) that \Lz x — xq\ 
is no smaller than \x — x$\ — L. In this case, 



inf Vr > C _1 U\x — xo\ - L) 2 — nL 2 ) 
Thus, if we set r$ := AnL, then whenever \x — xo| > r o we have that (I: 



Xo\ 



and 



Xo\ 



nL 2 > jq\x — X()\ 2 . Combining these inequalities with the bound above, we have that 



v R (x)>C- 1 \x-x \ 2 , VxeO\B'(4 \x\<R- 



Arguing in a similar manner we can obtain an upper bound on \vr\ for all points in Dr, so that 

\vr(x)\ <C{1 + \x~x \ 2 ), Vx 6 Dr. 

for a universal constant C. In particular, since this estimate independent of R we have enough compact- 
ness to do a Cantor diagonalization argument as R — > +00. Therefore there is a sequence Rk — > +00 
such that VR k converges uniformly on compact subsets of Oi to some v : 0\ —¥ K, which is continuous 
and belongs to Hl oc {0\). Moreover, for some universal constants C,C\,C% all > 1 we have 



|AS| 
d n v 
v 
v 



< C x 
= 

> C~ x \x 

< C%\x — Xq 



2 



in Oi, 
on dO u 
mO\Bl a {x»), 
\nO\13lixo). 



Setting clq := max{l, (C\C) }, a\ := C2C 1 we obtain the lemma with v := C 1 



a 



The problem of building a barrier now reduces to defining it inside B^ r (xq). To do this, we first 
observe there is a small ball around xq with universal radius where a standard parabola can be used. 

Proposition 4.8. Under Assumption [2] there exists a universal constant po such that for any xo € O, 
the function 



v(x) 



satisfies d n v > on (DO) n B Po (xq)- 



Proof. We pick r small enough so that for any x 6 0(u>) the portion of the hypersurface dO{w) lying 
in B r (x) is a radial Lipschitz graph with respect to x. That such an r exists follows from Assumption 
[2] and its size is determined by the constants in this assumption, in particular, r can be chosen as a 
universal constant. 

In this case, we see that for any x' £ B r (x) U dO(ui) x' — x is pointing towards the exterior of O at 
x', that is, d n v > on B r (x) n dO(w). □ 

If we had po > ro, where ro, po are as defined in the previous Lemma and Proposition, we could try 
to combing the function from Lemma |4.7| with the standard parabola. The question then is whether the 
resulting function is also a supersolution, that is, satisfies — Av > — 1 in O and d n v > on dO. However, 
we expect po to be much smaller than ro, so first we must extend the parabola from Proposition 4.8 to 
a barrier defined a ball of radius much larger than ro. 

Proposition 4.9. There are universal constants C > 1, Rq > tq and a e (0, 1) such that given e > 0, 
xo G O there is a semicontinuous function v G H 1 (B Roe (xo)) such that 



-Av 
d n v 
v 

V 



> -1 

> 



> \x- 



on dO, 
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where ro is as in Lemma |4.7| Furthermore, v satisfies the bound 

inf D-supu> C -1 , (4.1) 
The infimum being over B\ Rq (xo) \ Bl Ro / 2 ( x o) and the supremum over B\ ro {x ). 



Proof. As in the proof of Lemma [4. 7[ by scaling we only need to construct v for e = 1. Let Ro, a > 
be fixed and to be chosen later, and consider an obstacle problem in B Ra (x ) with obstacle given by 
— xq\ 2 and right hand side given by —2a. The solution, which we denote v, satisfies 

-AS >-2a mB Ro (x ), 
d n v > on dO, 

v >^\x-x \ 2 mB Ro {x Q ). 

Applying the theory from Section [3] leads to the bound 

sup v < C{^Rl + 2aRl) < CaR 2 ,. 

Indeed, this is guaranteed by the Stampacchia maximum principle (Lemma |3.14 ) and the fact that v 
must agree with the obstacle somewhere in B R {xq). In this case, pick a = a(C, Ro, po) such that 

CaR 2 < £(pg/2). 

Then, it follows that the function 



v(x) 



min {^k\ x ~ xo\ 2 ,v{x)} wB* o (x ), 



v(x) in B Ro (x Q ) \B x paJ2 (x ), 
is well defined. This is due to min{^-|a: — xo\ 2 , v(x)} — v(x) in Bp Q (xo) \ Bp g ^ 2 (xo), thanks to 

v(x) < ±( P o/2) 2 <&\x- x \ 2 in B^xo) \ B^xo). 
It is also clear that —Av > — 1 in B R (xq) and d n v > on (dO) n B Ro {xq). On the other hand, 

v(x) > ^\x-xo\ 2 VxeB Ro (xo). 



Finally, (4.1 1 can be obtain with an elementary barrier argument. □ 

Proof of Lemma \J-4\ We use again the scale invariance to reduce the construction to the case e = 1. We 
are going to "glue" the functions built previously by exploiting once again the fact that the minimum 
of two supersolutions is a supersolution. To this end, let j3 > to be chosen later, vi(x) the function 
constructed in Proposition 4.9 and v 2 {x) the function 

v 2 {x) := /3v 2 (x) + sup v\, 



where v 2 is the function constructed in Lemma 4.7 Observe that v± is a function defined in B Ro (xq) 
and v 2 is defined in O \ B^ (xo). Then, the barrier v is defined in O by 

!vi(x) inBl (x ), 
mm{v 1 (x),v 2 (x)} in B Ro (x ) \ B 1 ^, (4.2) 
v 2 (x) ^0\B Ro/2 . 

Let us see this v is indeed well defined. First, note that 

min-^i,^} = vi in Bl o (x ) \ Bl o/2 (x ), 
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which is guaranteed by the definition of v-i. On the other hand, due to (4.1 1 there is a /? determined by 
v\ , a,\ and R (and thus universal) such that 

(3a,iR 2 + sup Vi < inf v\. 

This guarantees that 

min{t;2,wi} = v 2 in B^(a;o) \ S]j /2 (x ), 

proving that u given by (4.2 1 is a well defined function. Furthermore, as v± and vi are supersolutions of 
the elliptic equation in their domains of definition, the same holds for v, thus 

-Aw > -1 in O, d n v > on dO. 

Next, note that V\ and v 2 satisfy the bound (in their respective domains) 

Vi(x) > a\x - x Q \ 2 , 

for some small universal constant a, thus v(x) > a\x — xq\ 2 in O. Finally, since vi(xq) — and v± = v 
near xq it follows that v(xq) = 0. This constructs the desired barrier when e = 1, finishing the proof. □ 

With the barrier at our disposal, we can now proceed to prove the non-degeneracy estimate. 



Proof of Theorem 4-3 We carry out the proof of the non-degeneracy estimate in two steps. Note that 
since the estimate will be uniform for x in the support of p c , we only need to prove it when 

xe{ P e (-,t) >0} 

Step 1. Suppose first that (x,t) and r are chosen such that 

r(-,t) > | mB*(x). (4.3) 



Let v be the barrier given by Lemma 4.4 for the point x and e > 0. The function w(-,t) — p e (-, t) — is 
strictly positive at x, which is guaranteed by the fact that p e (x, t) > and that v(x) = by construction. 
Then, for any r > and with D r :— B^(x) n {p e > 0} we have 

supio > 

D r 

Moreover, since Av < 1 in O e and d n v > on dO f l 

Aw > in Bfa), d n w < on dO e . 
As a consequence of the maximum principle, we obtain 

sup w = sup w > 0. 

d"D r D T 

Since v > everywhere, the supremum on the left cannot be achieved on the portion of dD r lying within 
d UJ {p c (-,t) > 0}, since there we have w = — v < 0. We conclude that the supremum over d^D r must 
be achieved at a point in d u B^,(x). In other words, 

sup [p=(., *)-!«] >0. 



Then, using that v(x) > ar 2 outside B^.(x), we get sup p e (-,t) > \ar 



2 



Bf(x) 

We have proved the following: for any x, t and r for which ( |4.3[ ) holds we have the estimate 

sup p e (-,t) > \ar 2 . (4.4) 

B«(x) 



5iep ^. We are left to prove the theorem without assuming (4.3). For any k G N write tk — \k and 
also to = 0. Arguing by induction, we will show for any k G N that if t G and x, r are such 

that r < d(x,D ), then 

sup p e (-,t) > a4-( fc+2 V 2 (4.5) 
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The theorem follows from this since for C = log 2 (a) + 10 we have e _c, ° (t+1) < a4 (fc+2) V k, t> t k . 
To prove this first for k = 1, note that (4.3| holds for any combination of x,t,r for which t < | and 



r < Dq), since in this case we have B^.{x) C (Dq) c , so 



{ P '>q} ds>\, in B*(x) 



Then applying the estimate from Step 1, we prove the base case of the induction argument. Next, 
suppose we have proved (4.5) up to k — 1. Then take any t G (tk—i,tk) and consider the function 
h(x) :— p £ (x, t)—p £ (x, tk~i). Recall that p e (x, t) is increasing in t, therefore h > 0, moreover, p £ (-,t) > h. 
Assume first that B*/ 2 (x) lies entirely outside {p c (-,ife_i) > 0}, then h solves 

ft 



Ah = l = 



ds > \ in D r , 



where D r := ({p e (-,t) > 0} \ {p e (-,tk-i) > 0}) n B^ 2 (x). Then taking v as the barrier from Lemma 
for the point x, we have again that w — h ~ has a positive supremum in D r , and that by the 



4.4 



maximum principle it must be achieved on dD r , and arguing again as in Step 1 we conclude 

sup h> sup h> \a{r/2) 2 = \ar 2 > aA- {k+2) r 2 . 

B' r (x) B' r/3 {x) 

and thus we get the same lower bound for the supremum of p e (-,t) over B e .{x). Now, if it was not the 
case that B £ , 2 (x) c {p t {- 1 tk-\) > 0}, we would know there is a point G {p e (-,tk-i) > 0} such that 

B' r/3 (x ) C Bl{x) 

and by the inductive hypothesis we would have 



sup p € (-,t k ^i) > a4 



-((fc-l)+2) 



B e r/2 {x ) 



(r/2) 



a4 



Since p e (x,t) is nondecreasing in t we get sup p c (-,ifc-i) > sup p e (-,ifc_i), so we arrive at 

sup p £ (-,i fe _i) > a4T (k+2) r 2 . 
B%{x) 

With this the induction argument is finished and the theorem is proved. □ 

Remark 4.10. Observe that the non-dcgcncracy coefficient in Theorem |4.3| decays exponentially with 
time. While this is likely not an optimal estimate, it is enough for our current purposes. 

4.2. Estimates with respect to time. Let us now discuss the regularity of p e in the time variable. 
Since f e (x,t) is non-decreasing in time it follows that p e is non-decreasing as well. Moreover, in [28] . 
it was straightforward from the formulation of the problem that there is some universal constant C 
independent of e such that 

\p s {x,h)-p e {xM)\ 



sup 



< C 



(4.6) 



This estimate still holds in our setting, despite the appearance of the extra term J Q * I{j,e(., s )>o} ds. 

Proposition 4.11. Almost surely there is a C > such that (4.6 1 holds for all < e < 1. 

Proof. Given t\ and ti note that p £ (.,ii) each solve obstacle problems in O t with zero data at infinity, 
zero obstacle and right hand sides given by f e (-,ti). Now, note that 



1/^0-/^2)1 = 



Hp'>0} 



ds — 



Hp>>o} 



< t- 



then, using a standard barrier argument, (4.6 1 follows as done for instance in [B]. 



□ 
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We can combine this proposition with Lemma 4.3 to yield the following estimate in time. 

Lemma 4.12. The positivity set of p e , that is D t (p e ) := {x \ p e (x,t) > 0} is evolving continuously in 
time. Namely, P-a.s. there is a universal C > such that if s < t and e € (0, 1) then 

B s (f) C V t (f) C D s (f) + B c ^- S . 

Proof. The first inclusion is due to the fact that p e is non-decreasing in the time variable. To show the 
second inclusion, suppose x E r t (p £ ) is such that Bg(x) n D s (p e ) — 0. It follows that 

p e (y,s) = for all y e B e s (x) 

Then \p e (y,t)\ < C(t - s) for all y G Bg(x), since x <E T t (p e ) Lemma 4.3 says that p e (y,t) > C8 2 for 
some y € Bg(x) which gives S 2 < C(t — s) and proves the lemma. 

As a corollary, we get continuity of the droplet for t — 0, which will be useful in Section [6} 

Corollary 4.13. P-a.s. there is a C such that Dg C D t (p c ) C D„ + B c ^- t for any t and e e (0, 1). 

5. Convergence of the obstacle problems 



□ 



Due to the theory in Section [3j the functions p e are uniformly continuous in compact subsets of 
( xM and bounded in H^ oc (O e ). Thus, P-a.s. from any sequence ej, — > we may extract a uniformly 
converging subsequence. In this section we identify the limit of this sequence as the unique solution of 
a homogeneous free boundary problem. 

Theorem 5.1. Let p : R d x M+ — ?► R be the (unique) weak solution to 

-div(AVp) = /i (-Id§ + Jo I{p(., s )>o} dsj hi {p > 0} 
p = |Vp| = on 9{p > 0} 



(5.1) 



Then, if {p e } e are weak solutions of (2.7) with initial data Dq, we have p e — > p P-a.s. in L p norm for 
every p € [l,oo] (see Definition 2.7). Here A is the effective diffusivity of O (see Section [7j Definition 



7.2) and jj. £ [0, 1] is the volume density of O given in Definition 2.4 



Remark 5.2. The notion of weak solution for (5.1 1 is the same as that given for (2.7) in Definition 2.17 
when O = R d . The existence theory for (5.1) above problem is discussed in Appendix |a} Uniqueness 
for this problem will later follow (Corollary 6.4) by combining the uniqueness results for the Hele-S haw 
problem in E d (see [53]) together with the theory developed first in 
which is independent of this section. 



and extended in Theorem 



6.1 



First, let us carry out the extraction of a converging subsequence in detail. 



Proposition 5.3. Consider weak solutions p e as in Definition 2.17 with some fixed initial data Dq. 
Then, given a sequence ej, > 0, — > we have that for P-almost every u> there is a subsequence e' k and a 



continuous function p : 
we have 



with compact support such that for any x G O e (uS), and R,T > 
\p4(x,t,u)-p(x,t)\=0. 



lim sup 

fc ^°° Bfj(s)x[0,T] 



Proof. For each e and P-almost every u> let p e : It 
agrees with p e (w) in O e {u). That such a p e exists 



4.11 



x M + — > K be a Holder continuous function which 
-a.s. is a consequence of Lemma 4.2 and Proposition 
which also tell us that the Holder norm of p e , and thus that of p e , is bounded uniformly for e € (0, 1). 
Then, given any sequence et -> we only need to apply the Arzela-Ascoli Theorem and a Cantor 
diagonalization argument to get a Holder continuous function p(x, t) such that p tk — > p uniformly in 
compact subsets of M. d 



that is 



lim sup \p tk (x, t, lu) 

fc - > °° B R (x)x[0,T] 



p{x,t)\=0 

for any R,T > 0. Since p £k = p tk in B e R (x) x [0,T] c Br(x) x [0, T] by construction, we are done. 



□ 
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As outlined in Section |2.5[ we will use the free boundary estimates from Section [4] to show that 

/ e Io £ ^M/ inL 2 (M d ), 
where f(x, t) is given in terms of the putative limit p(x, t) by 



f(-,t) = -I D c + / I{ p (., s )>o} ds. 
Jo 

From here, stochastic homogenization theory (see Section [7|) will imply Theorem 5.1 The weak con- 
vergence stated above will be guaranteed by the next lemma. For convenience we let D t (p e ) := {x | 
p e (x,t) > 0} as in Section |4j 

Lemma 5.4. The following holds P-a.s. in uj: fix t > 0, and consider a sequence — > such that 
the functions p ek (•, w, t) converge uniformly to a function p(x). Then for every 5 > there exists k such 
that whenever k > ko we have that D t (p ek ) lies entirely in a 6 neighborhood of {p > 0} n O tk and that 
{p > 0} fl O ek lies in a S neighborhood of D t (p £fc ). 

Proof. For brevity we will omit u> throughout the proof. For each e we define 

p e :— max _min \x — y\, 

xeD t ( P ') y eD t {p)nO e 

all we want to show is that p tk — ¥ as k — > +oo, let us omit the k for now and simply write e. Note that 
the maximum and minimum above are always achieved, and we may assume without loss of generality 
that p e > 0. Therefore, we may pick a point x e € D t (p e ) such that 

p(;t)=0mB e pc (x e ). 



By Theorem 4.3 we know that as long as p t < ^d(x e ,Do), there is a point x t g Bp e (x e ) such that 

p e (x e7 t)>e' Co ^pl 

As a parenthesis, we note that the argument when ^(x^Dq) < p e is straightforward, since Dq C D t (p e ) 
for all t > 0, so we will ignore this case. Now, by the Holder estimate from Theorem |2.7[ we have 

P e (-,t) > e' Co{t+1) pl - Cr a in B*(x e ). 
Taking r = [^e^ Co{t+1) p 2 e ] 1/a we obtain 

p\-,t)> \e~ c ^ + V pi in B*{x e ). 
Since p = in B c Pi (x e ) this means that 

l e - c ^pl<\\f^t)-p\\ L ^ 0t) . 

By assumption the quantity on the right is going to zero along e/c as — > 0, so the inequality above 
says that p ek — > as well. A similar argument handles the the quantity 

p e := _max min \x — y\, 

x£D t ( P )no e v eD t (p«) 

so one can also show that p £fe — > as —> 0. We only need to note that p inherits the non-degeneracy 
property of Theorem |4.3| from p e due to the uniform convergence, so we can exchange the roles of p tk 
and p in the previous argument. 

This proves the first part of the lemma: all we need to do for any given 5 is to take ko be such that 
max{p £fc ,p efc } < S for all k > k - To deal with the convergence of the integral, note that {p(-,t) = 0} 
is a measurable set and thus d{p(-,t) = 0} has Lebesgue measure zero. Combining this with the fact 
that the domains D t (p ek ) are uniformly bounded sets (using that D is bounded and Corollary 4.13), we 
conclude that the volume of a neighborhood of d{p(-,t) = 0} can be made arbitrarily small. The first 
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part of the Lemma implies that for any S we can pick ko such that if N$ denotes a ^-neighborhood of 
d{p(-,t) — 0} then as long as k > ko we have 



L 



O e ,(aj)nN 5 



I V fc (-, -,*)><>} V)>0} I 



The quantity on the right is bounded from above by the volume of N$, which can be made arbitrary 
small, and we conclude that the limit of the integral is zero, as we wanted. □ 

Corollary 5.5. Under the same assumptions as the previous Lemma, we have P-a.s. in u> that 

dx = 0. 



lim 

k— ¥oo 



{p'k (.,w,t)>0} {P( )>0} I 



Integrating the previous convergence estimate with respect to time, we conclude that the functions 
/ Cfc converge in a strong enough norm to what we want. 

Corollary 5.6. Let i, w, e^, p £k and p be as in Lemma [5. 4| Then, for every t > fixed 



This convergence takes place in the L 2 norm, meaning that 



Hp(-,s)>0} 



ds. 



lim \\r*{x,u,t) - f(x,t)\\L*(o ..) =0. 
Proof. Since the convergence is uniform in space for every fixed t, Lcmma|5.4|and Corollary |5 . 5| say that 



lim 

k— >oo 



0«. 



IV* >o> ^ *) ~ Vo> *) I 2 dx = V 4 G (°« T )' 



(5.2) 



Integrating in i, we obtain the following estimate 

^ /o et * Jo IV* >o } s ) - I {p>o > (*> s ) i 2 dsdx 



< T fofo ck I V»><» (*>*)" \ p >oy(x,s)\ 2 dxds. 



Then from (5.2 1 and the dominated convergence theorem the corollary follows. 



□ 



With these lemmas at hand, we may now prove the convergence theorem. 
Proof of Theorem\5.1\ First of all, for any vanishing sequence of e's and P-a.s. we can find a function p 



and a subsequence e^- such that p £fc (-,u;, •) converges to p by applying Proposition 5.3 If we show this 
p is the unique weak solution to (5.1) then this limit would be independent of the original vanishing 
sequence of e's and u> € Q, meaning p e converges to p and the theorem would be proven. 
Recalling that f ek (x,t,uj) — — I(p*)c + J Q I{ P ' k (-, s )>o} ds, we claim that 

^O tk (x)I{p'k >0 }(x,t)f' ik (x,t) ->> nI{ p>Q y(x,t)f(x,t) inL 2 , (5.3) 
where / = — Id<= + Jq I{ p (-, s )>o} ds. Indeed, for any cj> £ L 2 (R d ) we want to bound 



h = 



l (/i>o) f(M)^I lP>01 /(M)^ 



dx 
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After rearranging the terms, we see that 

h< [ \i {pCk>a} -i {p>oi \\r k (x,t)t\dx+ [ \n{x,t)ci>- f{ x ,t)\K >oy 4>\dx 

< ll/ £t (^)llooP {pefc>0} -I {p>O }lll /2 |l0l| 2 + 11/%,*) - f(;t)h\\4>h. 



The first term goes to zero thanks to Corollary 5.5 and the bound ||/ £fc (., £)||oo < C(i) for all efc. At the 
same time, the second term goes to zero due to Corollary |5.6[ It follows that 



lim / I„ L «. „, f Cfc (x,t)ch dx — lim / I„ I, f (x, t)d> dx. 
On the other hand, Theorem |2 . 2 1 says that 



lim / I_ ip dx = (i / ip dx 

e— >G ' e * 



7.1 



and conclude that 



for ip = I { , t ) >0 }/(") t)4>i an d (5.3 1 follows. In this case, we can apply Theorem 

(AVp, V0) dx = J fi \I D % + J l{p(., a )>0} I{p(-,s)>O}0 ^ V^ei? 1 . 

where A is the effective diffusivity defined in Section [7] This shows that the limit t) is the unique 
weak solution to (5.1 1 and we are done. □ 



To finish, we remark that combining the uniqueness of the limit and Lemma |5.4| we have proved that 
the free boundaries d u {p e (-,t) > 0} converge uniformly to d{p(-,t) > 0} for every t > 0. 

Corollary 5.7. For P-almost every uj, t > and for every 6 > there exists eo such that for any 
e < to we have that D t (p e ) lies entirely in a 6 neighborhood of D t (p) fl O e and D t (p) n O e lies in a 8 
neighborhood of D t (p e ). 

Remark 5.8. In general one cannot show that D t is continuous in time with respect to the Hausdorff 
distance. Intuitively, what may occur is two fingers of ^>t(p) contact each other at time t = to along 
some contact set C having positive n — 1 dimensional Hausdorff measure. Then at the next moment C 
instantly disappears and we may have a discontinuity of D t (p) at t = to- 

6. HOMOGENIZATION 



We have now established the convergence of p e in the L°° norm (see Definition 2.7 1. Using this and 
the relationship between p e and u e , we will investigate the convergence of u e and D e . 

6.1. The relationship between if and p e . The first goal is to justify rigorously the relation 

p e (x,t,uj) — / u e (x, s, iv)ds. 



Since w plays no role in the proof of this relation it suffices to prove it for a deterministic domain O for 
which Assumption [2] holds, even for a deterministic O we shall use the notation. 

Theorem 6.1. Let O C R d be such that Assumption p] holds, and D C O such that d u D = <9"LV 



Let p : O x K + -> R be a weak solution to (|2.7|) and (2. 



For each t > 0, define «(•, t) as the unique weak solution to the boundary value problem (see Lemma 



6.2 for the definition of u for non-smooth d u {p(-,t) > 0}) 

-Aw = 1 in {p(-,t) > 0}nO, 

d n u = on dO, (6.1) 
u = on 9 w {p(-,t) > 0}. 
Then the following is true: 
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(a) u(x,t) is a viscosity solution of (2.3) (see Definition 2.12). Moreover u is continuous at t = 0, 
i.e. u(x,t) converges uniformly to u (x), corresponding to the initial phase D . 

(b) u is equal to the left hand side time derivative of p: 

u(x,t) = drp(x,t) := lim for t > Q 

The proof of this theorem is close to that of Theorem 3.1 in [25], but we present a detailed proof for 
completeness. A key idea in the proof is showing that if a supersolution touches the function u{x,t) 
from above, then its integral with respect to time touches p from above at the same point. This will in 
turn show that classical supersolutions starting above u stay above it for all times. 

First, we carry out the construction of the function u in detail. 



Lemma 6.2. Let p be as in the statement of Theorem 6.1 then define u(x,t) by 

u(x, t) := sup 11(2;). 

Here the supremum is taken (for each fixed t) over all lower semi-continuous functions v : O 
which there exists some s £ (0, t) such that in the viscosity sense we have 



for 



-Av > 1 in {p(-, s) > 0}, d n v = on dO and {v > 0} C {p(-, s) > 0}. 
Then, u(x,t) is a lower semi-continuous function solving in the viscosity sense 

-Au = 1 in {p(; t) > 0}, d n u = on dO. 
Moreover, we have the pointwise bound 

r x p{x,i) < u(x,t), (6.2) 

and the coincidence of the positivity sets 

{u > 0} = {p > 0} and d u {u > 0} = d u {p > 0}. 

Proof. That u(x,t) is lower semicontinuous follows because it is essentially the limit of an increasing 
sequence of lower semi-continuous functions, we omit the details. Since p{-, s)/s is continuous in x for 
every s > and it is a weak solution of (2.7), we conclude that it solves 

-A(p(-, s)/s) < 1 in {p(-, s) > 0} n O and d n u > on dO. 

Then, and application of the comparison principle says that 

p{x,s)/s < u{x,t) for all s < t. 



By the continuity of p in time (Proposition 4.11 ) and the semicontinuity of u it follows that 

p(x,t)/t < u(x,t) Vt > 0. 

Thus u(x, t) > whenever p(x, t) > 0. Again by lower semi-continuity, u = on d"{p > 0}, so 

{u > 0} = {p > 0}, d^{u > 0} = d^ip > 0}. 

Also note that, by definition, 

-Au(-,t) = 1 in {p(-,t) > 0} n O and u n = at regular points of dOC\{p{-,t) > 0}. 

Lastly, observe that — Au(-,t) < 1 in the interior of O, thus —Au*(-,t) < 1 in O and it* < on dO, 
which finishes the proof. □ 

Next, we prove the key lemma that will allow us to relate test functions for u with those for p. 
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Lemma 6.3. Let u be given by Lemma 6.2 and <\> a classical strict supersolution (see Definition 2.10 



in some cylinder E := £>*(x ) x (t — Tq, to), for some Xo,t and with To € (0,1). Suppose that u < 
in S , u < </> on {<j> > 0} and {u > 0} n <9p E CC {</> > 0}. 

There is a small enough r > (depending on u, (f>) such that if E := B^(xq) x (f — t, to) and 

p(x,t) := p(x, t) — p(x,to — r), 3>(£,t) := / (f>(x,s)ds V(x,t)&T,, 

J to — T 

then we have the strict inequality p < <I> in {<!> > 0} n E. 

Proof. First, let us compare the pointwise values of u and p in E. To this end we observe that 

{u > 0} = {p > 0} = {p > 0}. 
Moreover, we know that p(-,t) solves (for each fixed t 6 (to — r, to)) 

-Ap = (-I D g + Jo* I {p(-, S )>0} rfs) 1 {p(-,t)>0}\{p(-,t -r)>0} + (/ t *_ r I{ p (., s )> } d») hp(;t -r)>0} & 

<9„p = on <90. 

In particular, since r < r < 1 we obtain J Q I{ p (., s )>o} ds < r < 1 outside {p(-,to — T ) > 0}. On the 
other hand, {p(-,t) > 0} is monotone increasing in time, so Id c = 1 in {p(-,io ~~ r ) > 0} c . From all this 
it follows that ^ 

f-lDg +^ l{p(-,a)>0} !{p(-,t)>0}\{p(-,to-T)>0} < 0. 

Then, since {p > 0} = {p > 0} and {p(-,t) > 0} is monotone increasing with t we conclude that 



-Ap < / I{p(.,a)>o} ds %(-,*o-ro)>o} < (*-*o+T)I{p(. ) t)>o} inO, dp = on 90. 

Jto-r 

The above being not only in the H 1 sense but als o in the viscosity sense, due to the solution p(-, t) £ H 1 
(and thus p(-,t)) being continuous, by Theorem 2.7 In this case (t — to + r)~ 1 p(- : t) is a (viscosity) 
subsolution for the same elliptic equation solved by u in the domain {p(-,f) > 0} = {u(-,t) > 0}. 
Applying the comparison principle we obtain 

p(x, t)<(t-t + t)u(x, t) V(i,()£Ox (t - r, t ). 

Recall that u < (j> in {(f) > 0} and {u > 0} n dp Eo CC {<fi > 0}, then it follows there is a 5 > such that 

u + 5 < (j) on (6£E)n{</i > 0}. 

Therefore for all (x,t) on this same set we have p(x,t) < (t — to + r)((j>(x,t) — S). At the same time, 

(t — to + r)((f)(x, t) — S) = / <j>(x,s)ds+ I 4>(x,i) — <p(x,s) — S ds. 

Jto—T Jto—r 

so picking r small enough (depending on S and the continuity of <p) we have 



(t - t + r)(<p(x, t)-6) < / 4>(x, s) ds = $(x, t) in E n {$ > 0}. 

Jto-r 

Then, using that p = outside {<& > 0} we obtain 

p < $ on <9^E , p < $ on <9£E n {$ > 0}. (6.3) 

On the other hand, as (j) is a strict classical supersolution in Eo and E C Eo, it can be checked that (in 
the sense) 

-A<j> > (%(., t )>o})t - l{0(-,t)>o} in s > 9 «</> > on (SO) n E. 
Integrating this with respect to t from t — t to t, it follows that 

f -A$(-, t) > -I m . :t)>0 } + I { 0(. >to _ T) >o} + J i * n _ r I{0(., s )>o} rfs in E, 

1 <9„$ > on (SO) n E. 
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Also because is a supersolution, {(/>(■, t) > 0} is monotone increasing in t, thus {</) > 0} = {$ > 0} for 
t > to — t and the inequality above may be rewritten as 

-A$(-,t) > ^l{*(-,t -r)>0}<= - !{*(., s)>0} rfs )^ l{*(-,t)>0}- 

Finally, in E we have I{p>o} < I{$>o}i thus 

-Ap(-,t) < (t-t +T)I m .. t)>0} . 

On the other hand, 

^~I{*(-,t — r)>0}<= + l{*(-,a)>0} ds^j l{*(.,t)>0} > (* — *0 + T )I{*>0} ■ 

From all this we conclude that 

f -A$(.,f) > -Ap(-,t) in #(x ) x {*}, 

\ d n $(-,t) > d n p(-,t) on(dO)x{t}. 
Thus, by the comparison principle, 

SUp ($ -p)(;t) < Sup ($-p)(-,t)<0 

Bl{x a ) d-'B 1 r (x ) 



where we used (6.3| to get the strict inequality. This shows that $ < p in Eg and the lemma is proved. □ 



Proof of Theorem \6.1\ Let us show that u as defined by Lemma |6.2| is a subsolution in the sense of 
Definition |2.12[ the argument for the supersolution property is similar. If u were not a subsolution we 
could find some cylinder E and a classical strict supersolution <p(x,t) such that u < <f> on &£E and yet 
u > <f> somewhere in E. 

Let to be the supremum of the times t' such that u < <f> within E D {t < t'}. Then cf> touches u from 
above at some point (xq, to) € E. By using the strong maximum principle for the Poisson equation with 
Neumann data, we can easily obtain a contradiction if we have p(xo,to) > or (f>(xo,to) > 0. Thus we 
only need to rule out the case where x € 9"{p(-,io) > 0} n 9 w {^>(-,io) > 0}. 

By looking at the first time when u fails to lie strictly below </>, we obtain to > 0, Xo € O and small 
numbers r, r £ (0,1) such that u < <j) in E = B^( x q) x (*o — T 07^o) an d {u > 0} n E C {0 > 0}. 
Moreover, by a small perturbation argument on <f> we may assume that the free boundaries of u and 
stay separated on d p T,, so that u and (/) comply with the assumptions of Lemma 6.3 in a cylinder Eo 
centered at (xo,to)- 

Due to the Lemma we conclude then that p < $ in a neighborhood of (xo,to), but this is a contra- 
diction since xq € 9"{p(-,to) > 0} n 9"{<^(-,io) > 0}, which means that p(xo,to) = <!>(xo,to). From this 
contradiction it follows that (xo, to) cannot exist and thus that u can never go above cj>, proving that u 
is a subsolution. 

To prove u is a viscosity solution of (1.1) all is left is checking that u and u both satisfy the initial 
data. Due to Corollary |4.13| we have that 

D o c{p(. > t)>Q}cD o + B t i /a (0). 

As we are assuming dDo = 9Do, it follows that D t converges uniformly to Do in Hausdorff distance. 
This proves part (a) of the theorem, it remains to show (b). To this end let us define for t > h, 

-i .s p{x,t)-p(x,t-h) 
u h (x, t) ■= . 

Since is monotone decreasing with respect to h, there exists a function u~ such that u^(x, t) converges 
pointwise to u~(x,t) as h — > 0. Due to the inequality above we see that u~ < u. Let us show that also 
u~ > u to conclude that u~ = u. 

Consider any lower semicontinuous function <f> with support in s) > 0} for some s < t such that 

-A<f> < 1 in {p(-, s) > 0} n O, d n 4> = on DO. 
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Then, from the comparison principle we obtain Uh > <j> for any h > t — s. From the definition of u it 
follows then that u < u~ , and the theorem is proved. 

□ 



Recall that when dO is C 3 and <9Dq is C 1 we know that uniqueness holds for ( 2.3 1 . In light of Theorem 



6.1 as well as Theorem 2.5 we obtain uniqueness for (5.1) too. We recall that existence for this problem 
is discussed in Appendix |A| 

Corollary 6.4. Suppose 0, Do and u are as given in Theorem 16. 11 and suppose in addition that dO is 



C 3 and <9Do is C 1 , then the function p(x, t) := J Q u(x, s)ds is the unique weak solution of (5.1 ) 



6.2. Proof of Theorem 2.6, From now on, u e (x,t,w) will be the specific viscosity solution of ( |1.1| 



given in terms of p e for P-almost every u> via Theorem |6.1| In particular, we have 

u e (x,t,uj) := d^p e (x,t 1 uj). 
We recall the notion of upper and lower semi-continuous limits, 

u*(x,t, oj) := limsup u e (y,s,uj), 

(y,s)eO e ,e^(x,t),0 

u*(x,t,uj) := liminf u e (jj,s,u>). 
(y,s)eO e ,s->-(x,t),0 

Note that u* and are well-defined a.e. in R d x [0, oo) x f2, since Assumption [2] guarantees that 

R d = limsupO e P-a.s. 

These limits, known also as half-relaxed limits in the viscosity solution literature, are convenient for 
two reasons: first the functions u e are defined in different spaces (so we cannot talk about the usual 
pointwise limit in a;), and second, we must deal with the possible lack of regularity of w c , which may 
even be discontinuous at later times. As we have mentioned already, this is a problem even for a fixed 
e. The lack of regularity means we cannot simply arrive at uniform convergence by using a compactness 
argument. However, if we can show that u* = the uniform convergence follows. 

In any case, in working towards the proof of Theorem |2.6| we wish to investigate whether it* = u* 



and whether both are equal to the unique solution of (1.2). To this end, let p : R x M + — > R be the 
unique solut ion to (5.1) and let u : R d x R + —> R be the corresponding viscosity solution to (1.2) given 



by Theorem 6.1 In this case we have that u = d t p for a.e. x,t and 

J -div(AVu) = n in {p(-,i)>0}, 
\ u = on d{p(-,t) > 0}. 



We already know that p e — > p (Theorem 5.1 1 uniformly, this implies at once that 

u € -> u in L* oc (B £ R (x ) x [0,T]) P-a.s. 

for any Xq,R,T and any p 6 [l,+oo). This does give a homogenization result of u e to u in the LP 
topology, but it does not clarify whether u e converges to u pointwise. Moreover, using Corollary |5.7| 
together with d^{p e > 0} = d u {u t > 0} and d{p > 0} = d{u > 0} we obtain the uniform convergence 
of the free boundaries in the Hausdorff distance. 

All that is left to investigate is the uniform convergence of u e itself. Note that what we want to verify 
is, roughly speaking, whether taking the time derivative and taking the limit e for p t commute. 

The plan is to exploit the comparison principle for the limiting problem to relate u to u* and . We 
begin with the verification of the elliptic equations that u* and satisfy. 

Lemma 6.5. For every fixed t > and in the viscosity sense, we have 

-div (AVu*) < n in R d , (6.5) 
-div(AVw*) > fj, in {it*(-,i) > 0}. (6.6) 
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Proof. We only prove (6.6) explicitly, as (6.5) is proven with the same method. Suppose (6.6) does not 
hold in the viscosity sense, then there exists (x ,io) € {w* > 0}, to > 0, and a smooth function <f> such 
that u(-, to) — 4> has a maximum at xq while div(AV</>)(xo) > — /i. By adding c + S(x — xq) 2 to 4>{x) with 
c = m(xo, to) — 4>{x$) and with sufficiently small i5 > 0, we may assume that the maximum of to) — 4> 
is zero and strict. Now, since <p is smooth there is some small r > such that 

div(AV</>) > — fi in B r (xo) and u(-,to) < on dB r (xo). 

Next, consider the function if : B r (x ) — > K determined by 

—div(A\7ip) = — /i in B r (xo) and 93 = on dB r (xo). 

By the comparison principle, 93 > in B r (xo), and thus — 4> has an interior maximum in B r (xo). On 
the other hand, for each e > consider the function ip e : Bf,(xo) — > K determined by 

Atp e = —1 in ^(x ), (jC^ = on (9C e and (/? e = (j) on <92?^(iro). 

Then Theorem 7.1 (see Section[7| says that P-a.s. as e — >■ <y9 £ converges uniformly in compact subsets of 
B^,{xq) to (p. Therefore if e is small enough u e — ip e has an interior maximum in B^.(xo) while u e < (f> = ip f - 
011 dB^ixo), contradicting the strong comparison principle. This contradiction proves (6.6). 

Next, we recall that u e (x,t) > t~ 1 p e (x,t) and that 

lim|b e (-,t)-p(-,t)||^ (0e) =0, Vt >0, 

from where it follows that u*(x,t) > t~ 1 p(x,t). Thus 

{u > 0} C {u* > 0}. 

Then, applying the comparison principle for the regular Poisson equation, the above yields 

u < u . 

In light of this, all that is left to do to prove Theorem |2.6| is that 

u > u* and T(u) = T(u). 
This is the content of the next proposition. 



□ 



(6.7) 



(6.8) 



Proposition 6.6. [28l Proposition 3] u* is a viscosity subsolution of (1.2), in particular, u* < u. 
The proof is essentially finished, we just gather the facts we have proven so far. 



Proof. Due to Lemma 6.5, u* satisfies div(AVu*) > — fi in K d . Therefore we only need to check that u* 
satisfies the desired property for subsolutions on the free boundary d{u* > 0}. 

Suppose u* is not a subsolution of ( 1.2 ). Then there exists a point (xq, to) G d{u* > 0} and a classical 
supersolution cj) of (1.2) in a cylinder £ := B r (xo) x (to — r, to] which touches u* from above at (xo,to), 

u* < <j> in E, u*(x ,t ) = ct)(x ,t ). 

Clearly this contradicts the maximum principle whenever (so, to) is not on d{u* > 0} or d{4> > 0}. 
Thus we only treat in detail the case where d{u* > 0} and d{(f> > 0} touch at (x , t ). Note that, due to 
(6.7) we have u < u* < <f> in E, and we conclude u crosses 4> from below in E at (xo,to), contradicting 
the fact that u is a viscosity subsolution of (1.2). 



□ 



Proof of Theorem 2^6 We apply Corollary |5.7| and conclude that 

r(«*) = r(o = r(u) = r(«), 

D(u*) = D(u») = D(u) = D(u). 



(6.9) 



Then, thanks to (6.8) and Proposition 6.6 we obtain 

u < < u* < u. 

In particular, if u is continuous we have u — u which from the above estimate means that = u* , this 
finishes the proof of the theorem. □ 
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Finally, let us point out a case where it is known that u is continuous in time. 

Lemma 6.7. If D is star shaped with respect to a ball B r (0), then so is D t (u). Moreover, u is Holder 
continuous in both space and time. 

Proof. First let us show that D t (w) is star shaped with respect to B r (0). For this, for given S > and 
xq G B r (0) let us consider the perturbed function 

U(X, t) = (1 + 5y 2 u((l +5)(X- X ) + Xq, t). 



Then u also solves (1.2) with the initial data u(x,0) < u(x, 0), due to our assumption. Therefore 
Theorem 2.4 yields that u < u in M. d x [0, oo), and we conclude since 5 > is arbitrary. 

It follows that r t (u) is locally Lipschitz continuous at each time. Since u(-,t) solves a uniformly 
elliptic equation in its positive set with Dirichlct bou ndary condition on T t (u), it follows that u(-,t) is 



Holder continuous in D t (p). Moreover, due to Lemma 4.12 and D 4 (it) being star shaped it follows that 
u and T(u) is also Holder continuous over time. 

□ 

7. Linear Stochastic homogenization revisited 

In this section we discuss the homogenization of linear elliptic equations in random domains O e , 
which was used in Section [5] to identify the limit of p e . We continue to make extensive use the notation 
introduced in Section [2j 

Let O be a stationary ergodic domain and let O be such that O = {u> | G 0(oj)} as in Definition 
2.3 and consider a sequence of functions f e € L 2 (B £ ) such that P-a.s. we have Io e / e - / to some 



/ G L 2 (D). Here we will be studying the sequence {v e } t given by the Dirichlet problem 





= / e 


in DDO e , 


d n v e 


= 


on <90 e , 


v e 


= 


on d^D. 



We are interested in identifying a deterministic limit as e — > as done in the uniformly elliptic but 
stochastic case in [351 Section 3] and [34| . The case of a stationary ergodic domain with Neumann data 
has been previously considered by Jikov (see [26 ). In all cases, the main issue is dealing with corrector, 
once the corrector problem is solved, the method pioneered by Tartar in the periodic case [H Chapter 
1] provides the proof of Theorem |7.1| The review of the aforementioned results is given below. 

Theorem 7.1. Suppose that O satisfies Assumptions [T] and [2] and let u e and f e be as above. Then 
P-a.s. as e — > + we have the weak limit 

l a u c fiu in L 2 (D). 



Here /i = P(O) (as Definition 2.4) and u G Hq(D) is the unique solution of 

div(AVu) = / in D, u = on dD, (7.1) 
whee A is a positive definite matrix (the "effective diffusivity" ) that we will describe in detail below. 
Definition 7.1. The space of stationary potential vector fields, L 2 ot (£l) C L 2 (f2,K d ), is defined as 

L 2 pot (n) := {v G L 2 (n;R d ) | E[v] = 0, v{x,u) := v(T x u) is a potential vector field P-a.s.} 
Given a measurable subset O C £1 we will denote by Lp 0t (f2, O) the vector space 

L 2 pot (n, 0) = {ve L 2 (0-R d ) \3v'e L 2 pot (n) s.t. v = v' in O}. 
We will also work with a space of stationary ergodic scalar functions which give Hl oc functions, 
H 1 ^) := {0 G L 2 (rt) | 4,{x,uj) := <f>{T x u>) G ff£ c (R d ) P-a.s.}, 
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and the restriction of these functions to O, 

H^fi.O) := {(f) E L 2 {0) | 3 <j>' £ H l (Cl) s.t. = 0' in O}. 



We shall now use the extension operator from Proposition |3.2| together with the harmonic extension 
to build an extension operator the probability space. Although this requires assuming that the domain 
is Lipschitz and well connected (that is, Assumption [2]), it allows for a proof of homogenization in 
stationary ergodic domains that is conceptually closer to the proof of homogenization for a uniform 
elliptic operator in M. d . 

Lemma 7.2. If the domain O satisfies Assumption [2j then there exists an extension operator 

Ext :H l (n,0) ->H\Sl), 

which commutes with the ergodic action in the following sense: take any u g and x G R d , then 

if we apply Ext to the restriction of u o T x to O we have 

Ext(u) o T x = Ext(ti o T x ). 

Proof. We define Ext by first defining it in for every u, and then lifting the construction back 

to the probability space. The problem is doing this construction in some canonical way that guarantees 
it will commute with T Xl and thus that it can be lifted at all to functions in Q. 

It is not hard to see that the harmonic extension is well suited for this. Thus, given u € O) 
let u(-,uj) denote its realization as a function in 0{uS), and define Ext w ({f(-, uS)) as the uniq ue fu nction 
agreeing with u{x, u>) in 0(uj) such that for every k it minimizes (with Pfe(w) as Proposition 3.3 1 

|Vtf dx 

among all functions v G H 1 (Pk(ui)) that agree with u(x, uj) in Pfc(a;)\P/ c (w). The existence of Ext w (u(-, u>)) 
follows by Proposition |3.2| which is how we are using Assumption [2] It also implies the bound 

||Ext w (u(; w))||^i (Bp(0)) < C*||M(-,a;)||^ 1(e , (0)) , 

which holds for any r, e > with the universal constant C given by Proposition 3.2 



□ 

Lemma 7.3. If the domain O satisfies Assumption [2] then the space Lp Qt is complete. 

Proof. Consider a sequence {v' n } n of functions belonging to Lp 0t (fi) such that v' n — > v in the L 2 (0; M. d ) 
norm. Then, by definition, for each n there is a function u n G - H 1 (ri) such that Vu„ = v' n in O, and by 
Lemma 7.2 we know that if we define functions u' n :— Ext(u n ) G 'H 1 (il) then 

Wu'Junn) < C|KII«i(o,o) ^ C|l^llL2 ot (n,o)' 

It follows that the sequence {u' n } n is weakly compact in - H 1 (r2). Let u be a limit point of this sequence, 
then, the convergence of the v' n says that (in the L? norm) 

v' n -> Vit. 

We conclude that Vu = Vq in O, which shows that v G L 2 ot and proves the lemma. □ 



Lemma [773| leads to the following Theorem, which corresponds to the stochastic analogue of the "cell 
problem" in periodic homogenization with uniform ellipticity. 

Theorem 7.4. For any £ G K d there exists a sequence {w|}{ e> o}, u>| : M. d x £1 — > K such that 

(1) w| G H^ oc (R d ) P-a.s. 

(2) Vw|(x,a;) = Vu^(^,u;), and Vw| is a stationary ergodic vector field. 
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(3) For any e and P-a.s. we have 

j Aiu| = in ° £ ' 
\ d n w\ =0 ondO t . 

(4) Kmw\{x) = (x,£) in Lf oc P-a.s. 

Proof. We note that for each £ there is a unique ^ £ L 2 ot such that 

E[I o (£ + %<;)]=0 Vv£L 2 pot . (7.2) 

Indeed, this follows since Lp Qt is complete by Lemma |7.3l in which case we may apply the Lax-Milgram 
theorem. In particular, note that the map £ — > is clearly linear. Moreover, by Lemma 7.3 there is a 
function xs,(x,w) (not necessarily stationary) such that 

Vx C (a:,w) = * c (r x a;) 

whenever cj e T^O, that is whenever x £ O{oj), and x^(-,w) £ Hl oc P-a.s.. Finally, we define 

w\{x,u) = x-£ + exi(f ,w) - a £ (w) 

where a e (w) is chosen so that 

/ w\(x,u) dx = P-a.s.V e > 0. 

It is clear that w| satisfies Properties (1) and (2). On the other hand, note that 

IoH(ViD 1 (i,w),t)(w))dP(w) =0 Vu £ Lp Qt . 

In particular, this holds for v — V</>, when <^> € - H 1 (f7). Letting (f>(x,uj) := <p(T x uj), and integrating the 
above identity with respect to x over a ball Br(0), we obtain 

/ / loH^fx.w), V0(x,w))dP(w) da; = 0. 
Equivalently (using Fubini's theorem), 

/ / {Vw 1 (x,uj),V4>(x,iu))dx d¥(u) = 0. 
Jn ^Bjj(o) 

Which holds for any </> obtained from a <j) £ - H 1 (i7). From here one can see that 

/ (Vw 1 (x,u}),'V<l)(x,uj))dx = 0, 

since i? > was arbitrary, this proves Property (3). To show Property (4) note that P-a.s. the sequence 
is bounded in H^ oc . Due to the ergodicity of V«;| it follows that 

Vio| ^E[e + V X |] = £ in L 2 (ft). 

This shows that iu| converges in Lf oc and weakly in H 1 to a function with gradient £. Since w| has 
average zero in B\ it also follows that this function must be precisely x ■ £, and the theorem is proved. □ 

Remark 7.5. It is helpful to think of w| as corrections to the affine functions x — > x ■ £ for £ € E d as to 
make them "harmonic" in the domain O e with zero Neumann data. They contain enough information 
about the problem to allow us to identify the limiting operator, this is the essence of the method 
developed first by Tartar to prove homogenization for a periodic and uniformly elliptic matrix A(x). 
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Definition 7.2. Observe that, given £ 6 M. d , the vector 

E [l 0c (x)Vwl(x)] 

does not depend on x or e > 0, and is in fact a linear function of £. This defines a matrix A by setting 

A£ :=E \lo e (x)Vwl(x)] (7.3) 

This A will be called the effective diffusivity. It is not hard to see that A is a symmetric matrix with 
nonnegative eigenvalues, we shall see that Assumption [2] guarantees A is indeed positive definite. 

Lemma 7.6. Suppose that O satisfies Assumptions [l] and [2j then there is a universal Ao > such that 



Proof. First of all, we note that Theorem 2.7 applies to the functions w|, and conclude that in any cube 
Qr := Qr n we nave w^(x) £ • x P-a.s. in the L°° norm, as in Definition 2.7 In particular, if we 



fix R = 1 then P-a.s. for all small enough e we have 

w\{x) < -0.5|f| in {a: e Qf | x-£< -1}, 
w|(x)> 0.5|f| in{x€ Ql | x-£> 1}. 

Now, due to Assumption [2j there exists a network of disjoint cubes within Q\ going from | x ■ £ < — 1 
to | x ■ £ > 1 (see Proposition IB. 5l). Integrating Vu>| • £ along each of these tubes we can find a constant 
c > independent of e and £ such that 



/ Vwl£dx>c \tf 



\ ii', • c, ' 'his 

'2! 

for all small enough e. Taking the limit as e — > 0, and using Theorem |2.2| we conclude that 

mIQiI(^,0>c |CI 2 , 

we conclude that all the eigenvalues of Ao are strictly positive, and the lemma is proved. □ 
The next lemma follows Tartar [4j Chapter 1]. It concerns the limit of (for any fixed <f> € C^°(M d )) 

o e VM e V</> dx 

as e — > 0. Note both I 0e and Vu e are at best converging weakly (separately) ase^O, so computing 
this limit is a priori non-trivial. However, by using the fact that Vu e is not just any vector field (it has 
a potential), we can get around the lack of strong convergence. This is where we use the "corrected 
solutions" u>|, in fact, the argument does not depart far from the periodic case from this point on. 



Lemma 7.7. Let A be as in Definition |7.2[ and suppose that e& is a subsequence going to zero and u tk 
converges in L 2 (D) to some u € Hq(D). Then, after passing to a subsequence (which we still call e^) 
we have that 

Iq. Vu e * AVu in L 2 {D) 
Proof. Let X denote the weak limit of l 0r Vu E ' . Given any cf> e C£°(R d ), £ € R d we will show that 



X ■ dx = / AVu • £0 dx 
which will imply that X = AVu. Let u>| be as above, we have 

div(I 0c Vu ek )wl k (j> dx = / Vu tk ■ (V<f))w\dx + I Vu €k ■ V '{wf)(j> dx 



On the other hand, since div(Io efc Vw^') = 0, we also have for any k 

= / u £k ■ Vw\ k dx+ <j>Vu ek ■ Vwt k dx 
JO,. Jo,. 
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Therefore 



-/ divflo \7u ek )wl k 4> dx = \ Vu tk ■ {\7(j))wl k dx ~ \ u tk V(j) -Vwf dx 
Then, due to the convergence of the different terms, taking the limit in the last identity gives 



div(X)(x ■ £)<t> dx = 
Integrating by parts on the left we arrive at 



X ■ (V(j>)(x-Q dx 



uVcj) ■ A£ dx. 



div(X)(x ■£)(/> dx = - [ div(X)(x ■£)((> dx - [ X-^dx-j uV<j>-A£, 

JR d JR d JR d 



dx. 



It follows that 



/ 



X-tfdx 



uV<j> ■ A£ dx 



and integrating by parts on the right and using the symmetry of A we finish the proof. 



□ 



Proof of Theorem 7.1 From the equation solved by u e we have that Io e Vw c is bounded in L 2 (R d ) 
independently of e. Moreover, due to Theorem 2.7 we know that along a subsequence the functions u e 
converge to a function u £ Hl oc (D). 
Then Lemma 



7.7 



says (after passing to another subsequence) that Iq c Vu €k — 1 AVu. Then, since 

o (t VM £i ■V<f)dx = J f e (j) dx 
we conclude after passing to the limit that 

A\7u -\7(j)dx = fcj) dx. 



Thus u is a weak solution of the homogenized Dirichlet problem. Since this problem has a unique solution 
we see that u does not depend on the original subsequence tk so u e converges to u. □ 

Appendix A. Existence of solutions to the auxiliary problem 



In this section we discuss the existence of solutions to the auxiliary problem (2.7). 



Theorem A.l. Given Dq C O e , there is a solution to (2.7| in the sense of Definition 2.17 



For the sake of simplicity we will omit the e superscript through this section. The one tool we will 
need has to do with solutions to the obstacle problem in stationary ergodic domains, which falls within 
the abstract theory of variational inequalities. See [33, Cha pter II] for a proof. 

Recall the set B e C M. d x introduced in Definition 2.5 as e is fixed in this section, we omit it. 

Proposition A. 2. Given any function / : B — » K there exists a unique function 

p[f] :B^R 

such that p[f] > 0, p[f] E H l (0) P-a.s. and 

V(p[/] - 4>) ■ V0 + [p[f] -4>)fdx>0 V ^ G H\0) P-a.s. 



Remark A. 3. We note that the property of p being a solution in the sense of Dcfinition |2.17| is equivalent 
to the identity 

p(;t)=p[.f(;t)} Vt>0, 

where f(-,i) = Id§ — / ^{p(-,s)>0} ds. Thus, the weak solution we are seeking can be thought of as a 
fixed point of the map above. 
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Proof of Theorem \A.1\ We construct the solution through an iterative procedure. Fix Do C O and 
define po by 

!Ap = -t in D , 
d n p = on dO(uj), 
po = inD^. 
For n > 1, define p n by p„(-, t) := p[f n (;t)] where /„ is given by 

/„(•,*) :=I D g - / I{ J) „_ 1 (., s )>o} ds. (A.l) 
Jo 

Let us show that f n +i < /n for all n. Arguing by induction, we show first f% < f\. Since pq > in D , 

h = I D . - t I Dn 

so that /2 < /i will be implied by p\ > a.e. in D . At the same time, this is equivalent to showing 
that {pi > 0} is dense in D (by the continuity of p\). If {p\ > 0} was not dense in D , there would be 
a small ball contained in Do where p\ = 0. Then, let <fio solve the Dirichlet problem in that ball with 
right hand side equal to — t, zero Neumann data on dO and zero Dirichlet data on the boundary of the 
ball. Then, integration by parts yields that (we extend (f>o by zero) 

V(pi - 4>) ■ + (pi - 4>)h dx<0 for <j) = Pl + O , 

which is impossible due to the definition of p\ . Thus p\ > a.e. in Do and f% < /i , as we wanted. 
For the inductive case, suppose that f n +i < fn for some n. Then, the comparison principle says that 

Pn+l = p[/n+l] > P[fn] = Pn- 

In particular 

l{p„>0} < l{p„ +1 >0}, 

and by integrating with respect to time this leads to f n +2 < fn+it as we wanted. 

We have shown not only that f n is monotone decreasing in n, but that p n is monotone increasing in 
n. Moreover, p n has a bounded norm H 1 norm P-a.s., and /„ is uniformly bounded pointwise for each 
t, so they must converge pointwise a.e. in B x to functions p and /. 

On the other hand, P-a.s. the functions p n (x, w, t) are uniformly continuous in compact subsets of 
0(lo) x R_|_, so they also converge uniformly for each fixed u to their pointwise limit p(x, ui,t). In this 
case, one may argue as in Section [5] to show that 

I {p>o} = limI{ p „>o} P-a.s., 
the limit being in L 2 sense. This shows that the pointwise limit of f n is in fact 

/ = fog - / I{ P (., s )>o} ds, 
Jo 

and we conclude that p(-, t) — p[f(-, t)] for all t > 0, so p is a solution and the theorem is proved. 

□ 

Appendix B. Pointwise bounds for the Green's function in O 
Here we obtain pointwise estimates for the Green's function G e of O e (u>), that is, the solution of 

AG £ (-,y,w) =*„(•) inO e (w), 
d n G\-,y,io) =0 on50 e M, 

which is a function G e (-, : 0(w) x 0(uj) — > M. These estimates say that for P-almost every u the 
function G e (x, y, to) is pointwise comparable with \x — y\ 2 ~ d , the Green's function of R d . 

Although we do not need them elsewhere in this work, bounds of this type are of interest in percolation 
theory, as hinted at above (see also the discussion in [3] as well as in [5]). They are also crucial in the 



e 
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analysis of boundary behavior of harmonic functions in Lipschitz domains [36j . Proofs of these bounds 
following De Giorgi's or Moser's approaches are currently not available (even in the discrete case). This 
motivates us to review the proof of the bounds for the Green's function in the context of a perforated 
domain, using Assumption [2j 

The main theorem of this section is the following. 

Theorem B.l. There exists a universal constant C > such that 

C-^x - y\ 2 - d < G e (x,y,uj) < C\x ~ y\ 2 - d , Vi,!/eO e (w) P-a.s. 

when d > 3, and replacing \x — y\ 2 ~ d above with log \x — y\ when d = 2. 

The proof is along the lines of the Littman-Weinberger-Stampacchia theory |36| for fundamental 
solutions of elliptic operators in divergence form. Beyond the elliptic theory in Section [3] and a capacity 
comparison estimate (Lemma B.6 ) the proof follows the classical one. In particular, the notion of capacity 
[TOl Chapter 4] will play a central role in what follows. 

Definition B.l. We define the capacity of a set E C W 1 by 

Cap(£T) = inf jy \Vv\ 2 dx , where v e H 1 ^ 4 ) and E C {v > 

Moreover, for e > and w e we introduce the capacity of E C K d with respect to O e (w), 

Cap, U (E) = inf \ [ \Vv\ 2 dx , where v € H 1 ^^)) and E n O e (u) C {v > 1}° I 

Definition B.2. Given E and any of the two notions of capacity, we can find a unique u g H 1 that 
achieves the infimum above. In either case this u will be called the capacitary potential of E. 

Remark B.2. It is immediate from the definition of capacity that if A c B then Cap(^4) < Cap(B), 
likewise Cap eui (A) < Cap ew (B). Moreover, we always have Cap(A) > Ca,p eui (A). 

The connection between capacity and the fundamental solution is illustrated by the next propositions. 

Proposition B.3. [See [Ml Section 7]] Given e>0,weO, a>0 and y eR d \ O e (ui) let 

J e (a,y,uj) := {x | G e (x, y, u) > a}. 

Then, 

C&p e ^(J e {a,y,uj)) = i. 
Proof. It is well known [19, Chapter 4] that the infimum is achieved by a function given by 

«(•)= / G £ (x,-, W )^(x), 



o,- 



where v is the so called capacitary distribution of the set J e {a,y,ui). In particular, one can check that 
v is continuous at y and that v(y) = 1. From this it follows that 



1= / G <L (x,y 1 u)dv(x). 
Jo, 

On the other hand dv is by definition supported on <9 W J € (a, y, u), and there we know that G e (x,y,uj) is 
identically equal to a , from the definition of J e (a,y,uj) (note we are using that G e is continuous away 
from y, which follows in particular from Section [3]). We conclude that 

1 = a\v\ = aCap £A ,(J e (a,y,w)) 

dividing both sides by a the proposition follows. □ 
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Proposition B.4. For any y £ K d \ O e (ui) and any r > we have 

inf G e {x,y,uj) < [Cap (B r (t/))1 < sup G e (x,y,uj) 

x£dB r (y) xedB r (y) 

Proof. We fix e, lu and y for the rest of the proof. Given r > 0, define 

a(r) = inf G e (x,y,uj), b(r) = sup G e (x,y 7 uj) (B.l) 

xedB r (y) xedB r (y) 

Then, it is clear that 

J e (b(r),t,u>) C B r (y) C J e (a(r),y,u>) 

and by the previous remark, 

Cap £iW (J £ (6(r),y,cj)) < Cap e ^(S r (y)) < Cap e ^( J e (a(r), y, w)) 
we then apply Proposition |B.3| and conclude the proof. □ 



The last part of the proof is modeled on a well known argument in percolation theory regarding the 
connectivity of the infinite percolation cluster, see for instance [27] or [53]. In our case, Assumption [2] 
guarantees that O is always "connected enough" so the proof is much simpler. We start by quantifying 
this connectivity. 

Proposition B.5. Fix L > l,e > and a point x$. Then P-a.s. for any x £ 8 u B\{xq) there is a 
Lipschitz path j(t) C O e such that for a universal constant C and e := (x — x )/\x — x \ we have 

inf dh(t),dO e ) > ed /3, 

t£[0,l] 

7(0)=*, 7(1) e ff-BiOro), 

| 7 (t)-(7(*),e)e|<Ce, 

/ |7(*)|<ft<CL. 
Jo 

Here do is the constant from Assumption [5] Moreover, for each e one can find almost surely a finite 
family of paths 7i,...,7jv satisfying all the properties above and such that d(7i)7j) > C*e for i =/= j. 
Moreover, if Tx, ...,T„ denote their corresponding Ce-tubular neighborhoods then 



Ti C O e (w) and 



> c 



Proof. Let cj be taken from the set of full measure for which 0(ui) satisfies the structure bounds, as 
explained in Assumption [2] Recall that in this case Pk(oj) denotes the separated components of 0{w). 
We start by taking the segment starting from x £ d u B\ in the direction e until we are edo/3 away from 
the closest obstruction Pfc(w). Consider the projection of this H e Pk(uj) on the hyperplane perpendicular 
to e. Then H e Pk(uj) CC B' e (0) and we can pick some r £ B e (uj) \ H e Pk(cu). 

Next, we move along the projection of the segment moving along r onto the boundary of the obstruc- 
tion until we reach at a point where e is supporting, then we move straight up until we are within do/3 
of the next obstruction. 

In each case, because the perforations are uniformly Lipschitz, after a uniformly bounded time we 
will have left the perforation, which shows we went up by a term of order e. Moreover, by picking C 
large enough, it is guaranteed that we will move again in the e direction before hitting the boundary of 
the cylinder, which guarantees the last inequality. Since we move by a uniform amount upward in finite 
time, it follows that we reach the top of cylinder in a controlled time. Lastly, the Lipschitz curve can be 
re-parametrized, so that t ranges form to 1. □ 
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Lemma B.6. There is a universal Co > such that P-a.s. given R > and x € O e (u>) we have 

Cap ejU (B£(a!o)) > c Cap(B i i(x )) V e > 0. 

Proof. Fix e > and w G fi such that the bounds in Assumption [2] hold. First let us suppose that R = 1. 
Let u e be the capacitary potential of B{(xq), then 

C & p eu (Bi(x ))= [ \Vu'\ 2 dx=[ \Vu"\ 2 dx 

JO, Jo c \B\(x ) 

Then, the homogenization result from Theorem |7.1| and the uniform Holder estimate from Theorem |2.7| 
apply to u e . This, together with a standard compactness argument, shows there is a universal constant 
L > such that 

u e < \ <mVBi (x ), 
in particular, this constant does not depend on e, xq ot lo. 

Given this Lq, Proposition |B . 5| says that for each e > there is a finite family of Lipschitz continuous 
paths 7i,.--,7iV (N depending on e) such that if T 1 ,...,T/v denote their e-tubular neighborhoods then 
they are pairwise disjoint, and for some Cq > independent of e we have 

|7i|~L Vi, 

JV 

and \jT t 

i=l 

Moreover, each 7$ connects a point in 23f (xo) with a point in {u < |}. Thus u £ (7i(i)) is equal to 1 for 
t = and less than 1/2 for f = 1, which says that 

-1 

Vw e ( 7l (t)) • 7< (t) dt > 1/2. 



> Cq. 



By applying Cauchy-Schwartz, this leads to 

CL / \Vu e ( lt (t))\ 2 dt>L Q /2 
Jo 

We can repeat this argument for the volume integral over Tj, and conclude for each Ti that 

/ |Vu| 2 dx > C- l \Ti\ 
for some larger but still universal constant C > 0. Since the Tj are all disjoint we have 

|Vu| 2 dz > c- 1 ! (Jt,| > c- l \Bi {x )\ > c- 1 ^, 



L 



IQ c \B\(x ) 

by making the universal constant C a bit larger. Letting c = C _1 LQCap(_B 1 (0))~ 1 , we conclude that 

c aP e ,w(Bi(a;o)) > c Cap(-Bi(a; )) 

Since the constant Co is independent of e > 0, we can reduce the case R 7^ 1 to this one by rescaling 
w £ ,Sf(xo) and e (u;), and the lemma is proved. 

□ 

All that is left is to gather the previous lemmas and prove the main theorem of this section. 



Proof of Theorem \B.1\ Using Remark |B.2| and Lemma |B.6| it follows that 

CMBAy))- 1 < Cap^OB^y))- 1 < c^Cap^d,))" 1 , 

which becomes 

-^^Cap^^^r^c--^. 
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Then, in light of Proposition |B.4| we conclude that 



inf G f '{x,y,u) < [Ca Pe JB e r (y))] 1 < sup G<(x,y,u). 

xed^B' r (y) - x£d"B' r (y) 



On the other hand, applying the Harnack inequality (Theorem 2.7 1 to G e (.,y) we obtain 



sup G e (x,y,uj) < C inf G e (x, y, ui), 

x£d"B' r (y) x£d"B' r (y) 

where C is universal. It follows that 

C- 1 -^ < inf <?(x )y> w) < sup G e ( W ) < Cc^-^. 

r xed"B*(y) x£d"B^.(y) r 

Equivalently, this shows there is some universal C such that for all e > we have 

I ° 1 2 <G e (x,y,^) < . P-a.s. 

\x — y\ a z \x — y\ a z 



□ 
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